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(N ■ ABSTRACT 

> ■ 

{Sj \ Determining temperatures in molecular clouds from ratios of CO rotational 

lines or from ratios of continuum emission in different wavelength bands suffers 
from reduced temperature sensitivity in the high-temperature limit. In theory, 
the ratio of far-IR, submillimeter, or millimeter continuum to that of a 13 C0 
\ (or C 18 0) rotational line can place reliable upper limits on the temperature of 

the dust and molecular gas. Consequently, far-infrared continuum data from the 
COBE/ DIRBE instrument and Nagoya 4-m 13 C0 J = 1 — > spectral line data 

^ | were used to plot 240 /^m/ 13 CO J = 1 — > intensity ratios against 140 /mi/240 fim 

dust color temperatures, allowing us to constrain the multiparsec-scale physical 

.^J | conditions in the Orion A and B molecular clouds. 

The best-fitting models to the Orion clouds consist of two components: a 
component near the surface of the clouds that is heated primarily by a very 
large-scale (i.e. ~ 1 kpc) interstellar radiation field and a component deeper 
within the clouds. The former has a fixed temperature and the latter has a 
range of temperatures that varies from one sightline to another. The models 
require a dust-gas temperature difference of 0±2 K and suggest that 40-50% of 
the Orion clouds are in the form of dust and gas with temperatures between 3 
and 10 K. These results have a number implications that are discussed in detail 
in later papers. These include stronger dust-gas thermal coupling and higher 
Galactic-scale molecular gas temperatures than are usually accepted, an improved 
explanation for the N(H2)/I(CO) conversion factor, and ruling out one dust grain 
alignment mechanism. 



Subject headings: ISM: molecules and dust — Orion 
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1. Introduction 



Development of a complete theory of star formation requires understanding the physical 
conditions of the molecular gas in which stars form, both on sub-parsec and multi-parsec 
scales. On the latter size scales, observations of molecular clouds can provide insights into 
large-scale patterns of star formation such as in spiral arms, rings, or in starburst regions 
of external galaxies (e.g., see Rand et al. 1992; Vogel et al. 1988; Aalto et al. 1999; Kuno 
et al. 1995; Wong & Blitz 2000; Sheth et al. 2000; Seaquist et al. 1998; Kikumoto et al. 
1998; Paglione et al. 1997; Giisten et al. 1993; Rieu et al. 1992; Wild et al. 1992; Harris 
et al. 1991; Tilanus et al. 1991; Turner et al. 1991; Verter and Rickard 1989; Lo et al. 
1987; Sorai et al. 2000; Garcfa-Burrillo et al. 2000; Harrison et al. 1999; Wall et al. 1991; 
Sage et al. 1991; Scoville et al. 1985; Nishiyama and Nakai 2001; Nishiyama et al. 2001; 
Kennicutt 1989; Martin and Kennicutt 2001; Rownd & Young 1999; Bryant and Scoville 
1999; Mulchaey and Regan 1997; Aalto et al. 1995; Wall et al. 1993; Puxley et al. 1990; 
Verter 1987, 1988). The diagnostics of molecular cloud physical, or excitation, conditions 
are usually obtained from observed molecular spectral lines. Despite being composed almost 
entirely of molecular hydrogen, only very warm few x 100 K) gas in interstellar molecular 
clouds is sampled in emission with the purely rotational lines of H2 (Rodriguez- Fernandez et 
al. 2001; Rosenthal et al. 2000; Wright et al. 1999; Parmar et al. 1991) and the ro-vibrational 
lines sample even warmer (>2000K) or fmorescently excited H 2 (e.g., Takami et al. 2000; 
Fuente et al. 1999; Moorwood and Oliva 1990, 1988; Black and van Dishoeck 1987). This 
is because the purely rotational transitions are quadrupolar and only occur between widely 
spaced (hz//k > 510 K) levels (for a more detailed discussion, see Wall 1999, and references 
therein). Consequently, estimating the properties of the bulk of the molecular gas, which 
is cold (kinetic temperatures of ~ 5-20 K, see, e.g., Sanders et al. 1985), requires observing 
lines of trace molecules with dipole transitions and closely spaced levels. The molecule CO 
and its isotopologues, like 13 CO and C 18 0, are the usual molecules of choice for estimating 
the large-scale properties of molecular clouds, because the most common isotopologue of 
CO - - 12 C 16 - - is the most abundant molecule after H 2 , the levels are closely spaced 
(hz//k = 5.5 K for J = 1 — > of 12 C 16 0), and the critical density of the lowest transition 
is low (i.e. n(H 2 ) ~ 3 x 10 3 cm -3 for J = 1 — > 0). These characteristics of CO ensure 
that the lines are often strong — with radiation temperatures or T R (i.e., Rayleigh- Jeans 
brightness temperatures as opposed to Planck brightness temperatures) of a few Kelvins 
- for 12 CO J = 1 — > observations of clouds in our Galaxy, (e.g., Sanders et al. 1985) 
and occur throughout most of the molecular ISM. Observations of lines of other molecules 
are also important for coming to a holistic picture of molecular clouds. For example, the 
rotational lines of CS sample the dense gas (n(H 2 )> 10 5 cm~ 3 ) most directly associated with 
star formation (Lada 1992). Nevertheless, because the rotational lines of CO sample the 
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excitation conditions in the bulk of the molecular gas, these lines provide insights into the 
environment in which the star-forming dense gas exists and, consequently, providing insights 
into star formation on the scales of many parsecs. Multi-parsec-scale observations of the 
lines of CO and its isotopologues, 13 CO and C 18 0, are carried out on external galaxies (e.g., 
see the references listed near the beginning of this paragraph) or by low-resolution (few 
arcminutes), wide-coverage (many degrees) mapping of Galactic molecular clouds (e.g., see 
Lee at al. 2001; Dame et al. 2001, 1987; Dame and Thaddeus 1994; Grabelsky et al. 1988; 
Bronfman et al. 1988; Cohen et al. 1986; Sakamoto et al. 1995; Robinson et al. 1988; Sanders 
et al. 1984; Tachihara et al. 2000; Mizuno et al. 1998; Kawamura et al. 1998; Heyer et al. 
1998; Boulanger et al. 1998; Onishi et al. 1996; Oka et al. 1998; Nagahama et al. 1998; 
Sakamoto et al. 1997a, 1994; Bally et al. 1987; Maddalena et al. 1986). 

Besides the usually observed main isotopologue 12 CO (written this way to explicitly 
distinguish it from the rarer isotopologues), the lines of the rarer isotopologues 13 CO (short 
for 13 C 16 0) and C 18 (short for 12 C 18 0) furnish very important additional information on 
molecular cloud structure. The lower-J lines of 12 CO are often optically thick and cannot 
probe deeply into the substructures, such as clumps or filaments, of molecular clouds. Also, 
the opacity of the 12 CO lines reduces their sensitivity to the effects of changing density and 
changing radiative trapping. Therefore, the lines of 13 CO or C 18 are crucial for realistically 
constraining the molecular cloud excitation conditions, because these lines, being from lower 
abundance isotopologues, are often optically thin on multi-parsec scales. This greater sen- 
sitivity of the 13 CO and C 18 lines to cloud structure and physical properties makes them 
capable of identifying interesting regions within molecular clouds. Of particular interest is 
locating regions of warm (e.g., ^ 50-100 K) molecular gas, because energetic events and mas- 
sive stars associated with recent star formation heat the surrounding molecular gas. Such 
gas is sometimes identified from ratios of strengths of CO rotational lines (e.g., Wilson et al. 
2001; Plume et al. 2000; Howe et al. 1993; Graf et al. 1993, 1990; Boreiko and Betz 1989; 
Fixsen et al. 1999; Harris et al. 1985; Harrison et al. 1999; Wall et al. 1991; Glisten et al. 
1993; Wild et al. 1992; Harris et al. 1991). 

Examples of how two of the line ratios involving the three lowest rotational lines of 
CO vary with gas kinetic temperature, T K , in local thermodynamic equilibrium or LTE are 
shown in Figure 1. All the curves in Figure 1 show the same trend: the line ratios lose their 
sensitivity to T K in the high-T K , or Rayleigh- Jeans, limit. This is a consequence of the line 
ratios depending on ratios of two Planck functions at different frequencies. The maximum T K 
to which a given line ratio is sensitive can be quantified by adopting a calibration uncertainty 
for the line ratios. If we adopt the optimistic value of 20% for the uncertainties in the ratios, 
then a ratio that has an observed value ~ 80% of its asymptotic value in the Rayleigh- Jeans 
limit is observationally indistinguishable from this asymptotic value. Similarly, the value of 



-4- 



the kinetic temperature at this 80% level is the maximum kinetic temperature for which the 
given ratio is sensitive. Two important conclusions can be drawn from the plots of Figure 1: 

1. Optically thin lines, like those of 13 CO, have a greater range of sensitivity to T K than 
do optically thick lines, like those of 12 CO. Specifically, T R (J = 2 — > 1)/T R (J = 1 — > 0) 
loses T K sensitivity above 46 K in the optically thin case and above 10 K in the optically 
thick case. 

2. Ratios involving higher- J lines have a greater range of temperature sensitivity than 
those with lower-J lines. For T R (J = 3 — > 2)/T R (J = 2 — > 1) the maximum T K values 
are 73 K and 13 K for the optically thin and thick cases respectively. These values are 
higher than the corresponding values for the T R (J = 2 — > 1)/T R (J = 1 — > 0) ratio, 
especially for the optically thin case. 

It is clear then that unambiguous identification of warm molecular gas requires observations 
of high-J lines of 12 CO and, sometimes, the optically thin lines of 13 CO (see references cited 
at the end of the previous paragraph). (It should be mentioned that simply observing the 
optically thick lines of 12 CO does not automatically give the kinetic temperature, because 
the observed lines are not necessarily thermalized and the gas does not necessarily fill the 
beam in the velocity interval about the line peak.) 

However, identifying warm gas in this manner can be complicated. One complication is 
that observing higher and higher lines up the rotational ladder becomes increasingly difficult 
with ground-based observations, because the atmosphere becomes increasingly opaque, on 
average, and because sub-millimeter receivers become increasingly noisier with decreasing 
observational wavelength. Another complication is that non-LTE effects become more im- 
portant higher in the rotational ladder: distinguishing between the density effects on the line 
ratios and kinetic temperature effects becomes difficult. One way to separate such effects is 
to observe many CO rotational lines, but this can be observationally intensive — especially 
if high-J 13 CO lines are also observed. Another approach to finding warm gas is using the 
12 CO/ 13 CO intensity ratio. This ratio is very sensitive to the molecular gas kinetic tempera- 
ture in the LTE, high-T K limit. For example, the 12 CO J = 1 — > 0/ 13 CO J = 1 — > intensity 
ratio in LTE goes like T 2 to within 25% for T K > 20 K. While the advantage of this ratio is 
that it only requires two lines with very similar frequencies, the considerable disadvantage 
is that the 12 CO and 13 CO lines occur in very different optical depth regimes. This results 
in three problems. One is that the 12 CO line will be effectively probing a different volume of 
gas from that of the 13 CO line. Another is that the levels of radiative trapping will be very 
different between the 12 CO and 13 CO lines, resulting in some sensitivity to density for the 
12 CO/ 13 CO intensity ratio. The third problem caused by the very different opacities in the 
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12 C0 and 13 C0 lines is that the 12 CO/ 13 CO intensity ratio will have some dependence on gas 
column density; indeed, in LTE this ratio is the inverse of the optical depth of the 13 CO line 
(for the 12 CO line in the optically thick limit and 13 CO in the optically thin limit). While the 
problems above are not necessarily prohibitively difficult, it would be helpful nonetheless to 
have two optically thin tracers with very different dependences on temperature, with relative 
insensitivity to gas density, and that are ubiquitous in the molecular ISM. The use of such 
tracers would complement the use of lines of CO and of other molecules in diagnosing the 
presence of warm molecular gas simply and unambiguously. 

Another way to trace molecular gas, and interstellar gas in general, is to observe the 
continuum emission from the dust grains found throughout most of the ISM. This contin- 
uum emission is usually observed at mid-IR, far-IR, submillimeter, and millimeter wave- 
lengths. On parsec scales, dust continuum emission is usually optically thin for wavelengths 
> 15 /im. All-sky surveys like IRAS and COBE/DIRBE cover specific IR wavelengths from 
the mid-IR to far-IR and near-IR to far-IR, respectively (IRAS Explanatory Supplement 
1988; COBE/DIRBE Explanatory Supplement 1998). Surveys at such wavelengths in the 
continuum are valuable for probing the structure and excitation of the ISM (see, for example, 
Dupac et al. 2000; Wall et al. 1996; Bally et al. 1991; Zhang et al. 1989; Werner et al. 1976; 
Heiles et al. 2000; Reach et al. 1998; Boulanger et al. 1998; Lagache et al. 1998; Goldsmith 
et al. 1997; Sodroski et al. 1994; Boulanger et al. 1990; Sellgren et al. 1990; Scoville & Good 
1989; Sodroski et al. 1989; Leisawitz & Hauser 1988). Ratios of continuum intensities at 
different wavelengths, for instance, can constrain estimates of dust grain temperatures. A 
physical model that includes dust-gas thermal coupling could then, in theory, constrain gas 
temperatures (e.g., Weingartner & Draine 2001; Mochizuki & Nakagawa 2000; Hollenbach 
1989; Tielens & Hollenbach 1985; Burke & Hollenbach 1983). Optically thin dust contin- 
uum intensity is T u B u (Td), where r„ is the optical depth at frequency z/, and B^(Td) is the 
Planck function evaluated at frequency v and dust temperature, Tj. In the far-infrared, the 
optical depth is often approximated as a power-law t u oc u 13 , where the emissivity index (3 
depends on the size and optical properties of the grains. Consequently, as with CO rota- 
tional lines, the ratio of continuum intensities at different wavelengths loses sensitivity to 
dust temperature in the Rayleigh- Jeans limit. As an example, we will consider the ratio of 
continuum intensities (per unit frequency) at the wavelength of 450 fim to that at 850 /im, 
i.e. Ij,(450 /im)/I^(850 /im), and adopt a 20% calibration uncertainty for this ratio. We will 
also adopt f3 — 2, which is a reasonable number for wavelengths longer than about 100 /im 
(Andriesse 1974; Seki & Yamamoto 1980; Wall et al. 1996). These adopted values imply 
that the 1^(450 /im)/I Jy (850 / um) ratio loses sensitivity to dust temperature for T d > 37 K 
(again from adopting 80% of the asymptotic value in the Rayleigh- Jeans limit). This max- 
imum Td can be improved by roughly a factor of 2 by combining the 450 /im observations 
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with observations at wavelengths longer than 850 /iin, but real improvement over this 37 K 
limit requires observations at wavelengths even shorter than 450 fim. Observations at such 
short wavelengths must deal with increased atmospheric opacity and sometimes must be 
space-based (e.g., IRAS and COBE). 

A more direct ground-based approach is suggested by the work of Schloerb et al. (1987), 
who show that a direct comparison of optically thin CO line emission with submillimeter con- 
tinuum can provide temperature estimates of the dust and molecular gas in molecular clouds 
(see also Swartz et al. 1989). Schloerb et al. (1987) derive an expression for the millimeter 
continuum to integrated C 18 J = 2 — > 1 line strength ratio, L,(submm)/I(C 18 J = 2 — > 1), 
which goes like T K Td in the high-temperature limit in LTE. For similar gas and dust tem- 
peratures, this means that the I !/ (submm)/I(C 18 J = 2 — > 1) ratio and, similarly, the 
I^(submm)/I( 13 CO J = 1 — > 0) ratio are actually more sensitive to temperature as that 
temperature increases. This is in stark contrast to ratios of rotational lines of a given iso- 
topologue of CO and to ratios of continuum intensities at different frequencies, which, as 
discussed above, lose sensitivity to temperature in the high-temperature limit. In addition, 
comparing the continuum intensity with that of a low-J line like the J = 1 — > line of 13 CO 
(or of C 18 0) reduces the sensitivity of the continuum-to-line ratio to gas density. In short, 
the intensity ratio of the submillimeter continuum to an optically thin CO line may be the 
sought-for alternative diagnostic of molecular gas temperature that was mentioned above. 

Despite the potential diagnostic advantages of the continuum-to-line ratio 
[e.g., L,(submm)/I( 13 CO J = 1 — > 0)], interpretation of this ratio and its spatial variations 
assumes the following: 

• The ratio of the mass of 13 CO to that of dust does not vary spatially. Even though the 
gas-to-dust mass ratio may vary, the 13 CO-to-dust mass ratio may not vary significantly 
within a single source, but may vary by up to an order of magnitude from source to 
source (Swartz et al. 1989). 

• The gas along a given line of sight is predominantly molecular. Appreciable quantities 
of atomic or ionized gas and their associated dust would complicate the interpretation 
of h,(submm)/I( 13 CO J = 1 — > 0). This is because the submillimeter continuum 
emission would not originate in dust only associated with molecular gas, but the 13 CO 
emission would originate only from the molecular gas. Consulting H I 21-cm and radio 
continuum data can check this assumption. 

• The molecular gas density does not vary spatially. If this assumption is false, then 
spatial gradients in the I^(submm)/I( 13 CO J = 1 — > 0) ratio may be probing spatial 
gradients in the molecular gas density rather than in the molecular gas temperature. 
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Again, observing the lower- J rotational lines of 13 CO, especially the J = 1 — > 0, line 
can minimize this problem. 

• The gas and dust temperatures have spatial gradients in the same direction. Even if 
T K 7^ T d , only assuming that the spatial gradients are in the same direction for both 
T K and T d is necessary. This is reasonable since both gas and dust are cooler with 
greater distance from heating sources, especially on multi-parsec scales. 

• Dust grain properties, like the dust mass absorption coefficient (k& v ) and emissivity 
index ((3), do not vary spatially. Observations of Orion, for example, show that such 
variations do occur for (3 on sub-parsec scales (e.g., Johnstone & Bally 1999), but, to 
first order, adopting a single value on multi-parsec scales is possible. 

• The 13 CO line used in the comparison is truly optically thin. If the 13 CO line opacity 
varies from optically thin to optically thick, the I^(submm)/I( 13 CO) ratio would again 
show gradients. The cause of the opacity change could still be due to a change in 
temperature, but it could also be due to a change in column density (and even volume 
density). Observing the corresponding C 18 line can test this. 

• There is no temperature gradient along the line of sight. Such a gradient would result 
in the submillimeter emission being biased toward the warmer regions along the line 
of sight and the 13 CO emission much less so biased (and even biased toward the colder 
regions if LTE applies). Even in the hypothetical case of T K = T d at every point 
along the line of sight, these biases could lead the observer, armed with observations at 
additional frequencies, to believe that T K ^ Td. This could also lead to the erroneous 
result that comparisons from one line of sight to another show no coupling between 
gas and dust temperatures. The effects of gradients along the line of sight can be 
minimized by observing clouds on multi-parsec scales, where such gradients would be 
less extreme. 

• Thermal and non-thermal gas emission make a negligible contribution to the submil- 
limeter continuum. This should indeed be the case. In M82, for example, Thronson 
et al. (1989) estimate that ~ 30% of the 1300 /xm flux in the central 30-45" is due to 
gas emission. This would be even less at shorter wavelengths and in regions with no 
dominant supernova remnants nor HI I regions. 

These assumptions may be reasonable first approximations for observations on multi-parsec 
scales or, depending on the assumption, can be tested for each source observed. 

Testing the overall validity of the above assumptions and, as a result, the behavior 
and effectiveness of the I^(submm)/I( 13 CO) ratio as a diagnostic of molecular cloud physical 
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conditions requires maps of a molecular cloud in the submillimeter or far-infrared (FIR) 
continuum, in the 13 CO J = 1 — > or C 18 J = 1 — > spectral lines, and of dust or 
gas temperatures, all with resolutions at multi-parsec scales. The most reliable maps of 
temperatures associated with molecular clouds would be from ratios of FIR maps. Ignoring 
certain complications for the moment, such as line-of-sight gradients in dust properties (e.g., 
Td), uncertainties in those properties (e.g., f3), and dust grains not in thermal equilibrium, 
the ratio of the intensities of FIR maps at two different wavelengths gives reliable dust 
temperatures because such intensities are in the Wien limit of the Planck function. In this 
limit, the dust temperature derived from the observed intensity ratio is insensitive to that 
ratio. For example, the observed continuum intensity ratio of 1^(140 fjm) j\ v (240 fjm) = 
1.66 would imply Td = 20 K. A 10% uncertainty in that intensity ratio would imply only 
a 5% uncertainty in the derived dust temperature (or ±1K). One complication, that of 
stochastically- heated dust grains (i.e., grains not in thermal equilibrium), can be effectively 
avoided by using two FIR maps at wavelengths of ^ 100 /im (for example, see Desert et 
al. 1990; Wall et al. 1996). This, therefore, rules out using IRAS maps, whose nomimal 
wavelengths do not go beyond 100 /im. On the other hand, maps from the DIRBE instrument 
aboard the COBE spacecraft cover the whole sky in ten IR bands from 1.25 //m out to 240 /im 
[COBE/ DIRBE Explanatory Supplement 1998). Two maps are at nominal wavelengths of 
longer than 100 /im (i.e., 140 /im and 240 /jm), making them virtually free of emission from 
stochastically-heated grains. Thus the DIRBE maps at 140 fim and 240 fim and their ratio 
can be compared with a map of 13 CO J = 1 — > (or C 18 J = 1 — > if available) of 
a given molecular cloud. Maps of entire clouds in our Galaxy in 13 CO J = 1 — > have 
recently become available (e.g., see Lee at al. 2001; Mizuno et al. 1998; Kawamura et al. 
1998; Nagahama et al. 1998). Therefore, the observational data exist to test the sensitivity 
and reliability to molecular cloud physical conditions of the ratio of the submillimeter or 
FIR continuum to 13 CO line intensities. 

The appropriate molecular cloud, or clouds, for which this testing should be applied 
must satisfy a number of criteria: 

- The test clouds should be bright in the desired tracers (i.e., FIR and 13 CO) and com- 
pletely or nearly completely mapped in these tracers, 

- largely mapped in tracers of atomic and ionized gas, such as HI 21-cm and radio 
continuum, 

- out of the Galactic plane to reduce confusion with foreground and background material, 

- several degrees in extent so as to accommodate many beams of the large DIRBE beam 
(0?7), 
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- and should have as large a range of dust temperatures as possible at the DIRBE 
resolution. 

The clouds of Orion most closely satisfy these criteria (see, e.g., Bally et al. 1991; Wall et 
al. 1996; Nagahama et al. 1998; Zhang & Green 1991; Green 1991; Chromey et al. 1989; 
Heiles & Habing 1974; Reich 1978; Berkhuijsen 1972; Haslam et al. 1970). While there 
are other clouds that may seem to satisfy these criteria as well (e.g., Chamaeleon-Musca, 
Mizuno et al. 1998), the Orion clouds have the largest range of dust temperatures at the 
DIRBE resolution (see Figure 2 of Lagache et al. 1998). The Orion A and B molecular 
clouds are bright in 13 CO J = 1 — > (Nagahama et al. 1998), and in the FIR continuum 
at 140 yum and 240 /zm (Wall et al. 1996), and has been completely or largely mapped in 
these tracers. Much of the Orion constellation has been mapped in HI (Zhang & Green 
1991; Green 1991; Chromey et al. 1989; Heiles & Habing 1974) and radio continuum (Reich 
1978; Berkhuijsen 1972; Haslam et al. 1970), thereby permitting tests of whether or not 
molecular gas dominates the cloud surface, or column, density for each line of sight. Much 
of Orion has also been mapped in 12 CO lines (Wilson et al. 2005; Lang et al. 2000; Ikeda et 
al. 1999; Sakamoto et al. 1997a, 1994; Maddalena et al. 1986), which can also provide crude 
estimates of the molecular gas column densities. Note that only Maddalena et al. (1986) 
have mapped the Orion all three major clouds: Orion A, Orion B and the A Orionis ring. 
The Orion clouds are about 15-20° out of the Galactic plane and about 6° x 3° (or about 
50 pc x 25 pc for an adopted distance of 450 pc) in size, down to some ill-defined zero level. 
An additional advantage to the clouds of Orion is that these clouds have been surveyed at 
many other wavelengths on the scales of many parsecs, thereby providing further checks on 
models. These other wavelengths include, for example, the optical surveys of the Orion OBI 
and A Ori OB associations (Brown et al. 1995, 1994; Warren & Hesser 1978, 1977; Murdin & 
Penston 1977; Blaauw 1964), the near-IR to far-IR all-sky survey of COBE/ DIRBE (Wall et 
al. 1996), the mid-IR to far-IR all-sky survey of IRAS (Bally et al. 1991), and submillimeter 
maps of small sections of Ori A and B (e.g., Johnstone et al. 2001; Mitchell et al. 2001; 
Johnstone & Bally 1999). Therefore, the Orion A and B molecular clouds were chosen for 
testing the behavior of the I l/ (submm)/I( 13 CO) ratio, using the DIRBE maps at 140 /im and 
240 pim and the 13 CO J = 1 -> data of Nagahama et al. (1998). 

In addition to testing the diagnostic abilities of the continuum-to-line ratio, the current 
data may also permit us to test theoretical models describing interactions between dust and 
gas. Specifically, models of dust-gas thermal coupling use very simple assumptions, such 
as large, single-size, spherical dust grains with smooth surfaces (e.g., see Goldsmith 2001; 
Burke & Hollenbach 1983). Are these simplistic assumptions sufficient for explaining the 
observations? This will be examined. 
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The current paper, or Paper I, describes the one- and two-component modeling of the 
data. Paper II (Wall 2006) in this series confirms and, in some cases, modifies the model re- 
sults using simulated data. Paper III (Wall 2006a) examines systematic effects not explicitly 
considered in the models and also discusses the scientific implications of the results. 

2. Data Processing 

2.1. DIRBE Data 

Details of the COBE mission and DIRBE instrument can be found in Boggess et al. 
(1992), Silverberg et al. (1993), Hauser et al. (1998), and the COBE/ DIRBE Explanatory 
Supplement (1998). Much in this section has been discussed in Wall et al. (1996) (hereafter 
referred to as W96) in Section 2.1 of that paper. The data in the current paper are from 
the 1998 release, whereas the data of W96 were of the 1993 release. Also, the current paper 
concentrates on Bands 9 (i.e. A = 140 /im) and 10 (A = 240 /zm), while W96 discussed all 
10 wavelength bands. Consequently, the data processing in this section largely concerns the 
maps in the 140 /im and 240 /im bands from the 1998-release of the DIRBE data. 

While a number of DIRBE data products are available, the ZSMA (Zodi-Subtracted 
Mission Averaged) maps were the most relevant for the current research. (See the COBE/ DIRBE 
Explanatory Supplement 1998, for more details of the DIRBE data products described in 
this paragraph.) These maps are the average sky brightnesses determined from averaging 
over the full 40 weeks of the COBE mission, thereby maximizing the signal-to-noise ratio 
of the data. Maximizing the signal-to-noise ratio is particularly important for the 140 /im 
and 240 /im bands, because these have the noisiest data. The average sky brightness in each 
wavelength band was determined after removal of the DIRBE interplanetary dust model for 
each of the ten wavelength bands, which effectively removes the zodiacal emission from the 
maps. The ZSMA maps are available in the form of FITS tables that include the pixel num- 
ber, the intensity in that pixel in MJy/sr, and the noise in that pixel in MJy/sr. The FITS 
tables were converted into maps in the equal-area Mollweide projection using IDL routines 
available from NASA/GSFC through the worldwide web. The fxbopen and fxbread routines 
were used to open the FITS files and read in the data, respectively. The data included the 
pixel numbers, sky brightnesses, and noise levels. The pix2xy routine was used to create 
maps from the data in the DIRBE sky-cube format. This was followed by converting the 
sky-cube maps to maps in the Mollweide projection using the reproj routine. Maps of both 
the sky brightnesses and of the rms noise levels were created this way. The maps of the rms 
noise levels are hereafter referred to as the sigma maps. 
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After the creation of the Mollweide maps using the IDL routines from NASA/GSFC, 
the section of the all-sky maps corresponding to the Orion region was extracted. As was 
done in W96, a cosecant-law background of the form a c csc(\b\) was subtracted from the 
sky-brightness maps, so as to remove foreground and background emission associated with 
the Galactic plane. The fitting procedure was virtually identical to that used in W96 with 
a slight change in the algorithm used to determine the Galactic latitude b from the pixel 
position. The earlier paper used a simplistic linear interpolation for determining b in the 
cosecant-law fitting routine. The current paper uses the exact expression that was used in 
creating the maps in the Mollweide projection. Nevertheless, over the latitude range of the 
Orion maps (i.e. b = — 5 to —31°), the improvement in the algorithm only contributes a 4% 
change to the scale factors, a c . As stated in W96, the uncertainties in the map intensities 
due to the uncertainty in the background subtraction is about 10%. The effects of this 
subtraction on the results will be discussed in Paper II (Wall 2006). 

These Orion maps were compared with those of the 1993 release of the DIRBE data. 
The brightnesses of the 140 and 240 /im maps were reduced by about 15% with respect to 
those of the 1993 maps made. Accordingly, luminosities and masses determined from these 
maps are also reduced by this amount. However, the 140 /xm/240 /im dust color temperatures 
are virtually unchanged. 

Before creating maps of the 140/im/240 /im dust color temperatures, the 140 /im and 
240 /im maps were smoothed, as was done in W96. Without smoothing these maps, the 
resultant dust temperature map appears to be very choppy: noticeable, large, and random 
variations that occur within a beam size. The maps were smoothed with a square flat-topped 
box that was 3 pixels x 3 pixels in width. Since the DIRBE beam size is roughly the same 
size (about 2.3 pixels x 2.3 pixels) and shape, the equivalent resolution after smoothing is 
roughly a/2 larger than the original 42' or about 1? The resultant maps of L,(140 /im) and 
L,(240 /xm) are shown in Figure 2. 

As with the intensity maps, the corresponding sigma maps must also be smoothed. 
Here, however, the smoothing procedure is not so straightforward. Suppose the smoothing 
box is n x pixels x n y pixels and Iij is the intensity at pixel ij in the unsmoothed map. Then 
the intensity in the smoothed map at the central pixel of the smoothing box, I s , is 



where N = the total number of pixels in the smoothing box. If the pixels are 

independent, then the rms noise level in the smoothed map at the central pixel of the 




n 
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smoothing box, a s , is 

*. 2 = ^££4- • ( 2 ) 

However, if the iy are already convolved with a beam the same size and shape as the 
smoothing box (i.e. n x pixels x n y pixels), then a s can be approximated by 

i n x riy 

1=1 J=l 

Equation (2) states that the noise in the smoothed map is a factor y/N below the quadrature 
average of the noise levels in the smoothing box. This is expected when the pixels are 
independent because the random walk applies. On the other hand, equation (3) states that 
the noise in the smoothed map is only a factor of \/2 below the quadrature average of the 
noise levels in the smoothing box. This is because the area of the effective beam is roughly 
doubled (exactly doubled in the case of Gaussian smoothing and a Gaussian beam) after 
convolution with a smoothing box the same size as the original beam. Equation (3) has been 
tested with simulated maps of random noise. After discarding pixels along the map edges, 
it was found that for 99.9% of the pixels equation (3) must be either corrected downward or 
upward by up to 30%. This is not a large correction to random noise estimates. Therefore, 
equation (3) was used in smoothing the sigma maps. 

The map of the dust color temperatures (which is the actual dust grain temperature for 
the large grains in thermal equilibrium in the case of homogeneous dust properties along the 
line of sight, assuming the spectral emissivity index is correct) is created in the same way as 
in W96. The dust temperature, T d , is derived from inverting 

n = ^ 9 B u9 (T d )K u9 (T d ,P) 

^ B uw (T d )K Ul0 (T d ,P) ' 

where TZ = l u9 /l Ul0 , l u9 and 1 V10 are the quoted DIRBE intensities at frequencies v 9 and 
u 1Q , v 9 is the frequency at A 9 = 140 jum, u 10 is the frequency at Ai = 240 /zm, B„(T d ) is 
the Planck function, and K u is the DIRBE band color correction. Even though these color 
corrections never vary by more than 8% from 1.0, their ratio, K v9 jK Vuv can vary by up to 
14% from 1.0 in the temperature range, Td = 15-30 K. Consequently, the color corrections 
were necessary for obtaining accurate dust color temperatures. Although there is evidence 
that (3 can vary spatially by 10-20% on sub-parsec scales (e.g., see Johnstone & Bally 1999), 
adopting (3 = 2.0 is still valid on the scale of a few parsecs for the reasons given in W96. 
Equation (4) was solved for T d by plotting T d as a function of 1Z over the temperature range 
2.1 K to 210 K, and fitting polynomials over three subranges within the full range. The 
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polynomial fits have a typical accuracy 0.08 K and never deviate by more than 0.4 K, which 
occurs at T d = 2.1 K and 45 K. The polynomial fits quickly determine T d for each pixel. The 
T d map is shown in Figure 3. 

The 1998 release of the DIRBE data includes the sigma maps and this allows the 
computation of the map of random errors in T d , i.e., a o"(T d ) map. To compute this map 
rapidly and easily, the Wien approximation is applied to the Planck functions in equation 
(4). The Wien version of (4) is differentiated and rearranged to give the following simple 
expression: 

<T d ) = T d o(^) 
T d A TZ ' 1 ' 

where it was assumed that the color corrections are only weak functions of T d and where 
A = h(z/ 9 — ^io)/k = 42.8 K. The constants h and k are Planck's constant and the Maxwell- 
Boltzmann constant, respectively. The relative error in 71, a(1Z)/1Z, was computed from 
the quadrature sum of the relative errors of the intensities in the 140 /im and 240 /im maps. 
It must be emphasized that even though the uncertainties in T d are computed from an 
expression that is based on the Wien approximation — i.e., (5) — no such approximation 
was used in determining the T d itself. Therefore, the T d that appears in (5) was computed 
from (4) (or an accurate approximation of it), which includes the full Planck function and the 
color corrections. Using the accurately determined T d improves the accuracy of the derived 
<r(T d ) values. A detailed comparison with cr(T d ) determined from a numerical differentiation 
of (4) shows that the cr(T d ) derived from (5) only needs to be corrected upwards by 20-25% in 
the temperature range 15-30 K. Such a correction was not applied and makes little difference 
to the results. (This is discussed later in more detail.) 

Besides the random uncertainties due to noise, there is the overall calibration or pho- 
tometric uncertainty of the map in each wavelength band. These are 10.6% and 11.6% 
in the 140 and 240 /im bands, respectively (Hauser et al. 1998). How these systematic 
uncertainties affect the derived color temperatures, T d , depends on the relative calibration 
uncertainty of the 140 /im band with respect to the 240 /im band. In W96 the calibration 
uncertainty in the 1^(140 /im)/L,(240 /jm) ratio was specified as 6%. Adopting that value for 
the current paper gives an uncertainty of ±0.4 K for T d = 15 K and ±1.6 K for T d = 30 K. 
These uncertainties are small compared to the relevant range of values of T d (i.e., 15-30 K) 
itself. Consequently, the systematic uncertainty in the 140 /im/240 /im color temperature is 
roughly equivalent to an upward or downward shift of about 1 K to the entire T d map. 
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2.2. 13 CO Data 

The 13 C0 data were obtained with the Nagoya 4-m antenna and cover both the Orion A 
and Orion B clouds. The beamsize is 2(7 and Orion A was mapped with a 2' grid spacing. Full 
details of the Orion A observations are given in Nagahama et al. (1998). The Orion B cloud 
was mapped with an 8' grid spacing. The data consist of velocity-integrated intensities of the 
13 CO J = 1 — > line and were gridded to the DIRBE resolution in the Mollweide projection 
to allow comparison with the DIRBE data. The regridding procedure was a straightforward 
conversion from the galactic longitude and latitude of the central position of each pointing of 
the beam on the sky to DIRBE pixel coordinates in the Mollweide projection. These pixel 
coordinates were rounded off to the nearest whole number to give the DIRBE pixel that 
contains the center of the beam. This simple approach does not take into account the cases 
where a given pointing of the 2'7 beam spills over the edges of given DIRBE pixel. However, 
this is not an important consideration given that the DIRBE-gridded map is then smoothed 
by 3 pixels x 3 pixels. Each DIRBE pixel was "hit" by a few pointings of the Nagoya beam 
and the corresponding intensities were averaged for each DIRBE pixel. The 1-a noise level 
of the Nagoya integrated intensity map was 0.7 K • km • s _1 . After regridding to the DIRBE 
resolution, the noise-level was computed for each DIRBE pixel by dividing the Nagoya map 
noise-level by the square-root of the number of hits of the Nagoya beam in that pixel. Hence, 
the regridding process produced both an integrated intensity map of the 13 CO J = 1 — > 
line and a corresponding sigma map. As was done with the DIRBE maps, the 13 CO map 
was smoothed by 3 pixels x 3 pixels. The sigma map was smoothed using equation (2), since 
the pixels are independent in this case. 

As mentioned above, the Orion B section of the original Nagoya map has an 8' grid 
spacing even though the beam is only 2'7 in size. This spacing gives gaps in the coverage. 
These gaps do not result in large systematic errors in the velocity-integrated intensities in the 
final map, largely because all the data were gridded to the DIRBE resolution. Systematic 
errors in the map intensities due to incomplete coverage can be estimated by comparing the 
Orion A observations on the grid with 2' spacing with those same observations on a grid with 
8' spacing. Both the 2' grid and 8'grid Orion A data were regridded to the DIRBE resolution 
and then smoothed by 3 pixels x 3 pixels. The median of the intensities of the map that was 
originally on an 8' grid were 10-11% lower than that of the map that was originally on a 2' 
grid. Therefore, it is possible that the Orion B observations need to be corrected upwards 
by, on average, about 10%. Such a correction was not applied, because this is only a small 
change that will not affect the results significantly and because application of a single scale 
factor to all the Orion B intensities is overly simplistic. 

A final consideration for comparison of the 13 CO data with the DIRBE data is whether 



- 15 - 



background subtraction, as was done for the DIRBE data, is necessary. As was found in W96 
for the 12 C0 data, such background subtraction for molecular line data produced virtually 
no change in their intensities in the Orion fields. Consequently, no background subtraction 
was applied to the 13 CO map. The map of velocity-integrated radiation temperature of the 
13 CO J = 1 — > line, or I( 13 CO), is shown in Figure 4. 



2.3. Other Datasets 

Other datasets used in this research include the map of the integrated intensities of the 
12 CO J = 1 — > line, the map of the integrated intensities of the H I 21-cm line, and the map 
of the peak 12 CO J = 1 — > radiation temperatures. The maps of the integrated intensities 
of the 12 CO J = 1 — > line and of the HI 21-cm line were described in W96 (see also 
Dame et al. 1987; Maddalena et al. 1986; Heiles & Habing 1974). As in W96, an integrated 
CO line intensity to molecular gas column density, or X-factor, of N(H 2 )/ / T H dv = 2.6 x 
10 20 cm -2 ■ (K • km ■ s -1 ) -1 is assumed. Recent work by Dame et al. (2001) suggests that 
an X-factor of ~ 2.2 x 10 20 cm -2 • (K • km ■ s -1 ) -1 is more appropriate for the Orion clouds. 
As will be discussed in Section 3.6, this is indeed a more appropriate value (at least for the 
one-component models fitted to the data), but the value of 2.6 x 10 20 cm -2 • (K • km • s^ 1 ) -1 
is adopted for consistency with W96. After regridding the data to the DIRBE resolution 
and performing the usual 3 pixel x 3 pixel smoothing, the noise levels in the CO and H I 
integrated intensity maps were estimated from the negative tails of the histograms. The 
1-cr noise levels are 0.7 and 0.2 in units of 10 20 cm~ 2 for the CO and HI, respectively. As 
described in W96, a smooth background representing large-scale background and foreground 
Galactic emission was subtracted from the HI map. No such background was subtracted 
from the CO map because it made little difference to the resulting intensities. Both the 
regridded CO and H I maps can be seen in Figure 3a of W96. 

The peak 12 CO J = 1 — > radiation temperatures were provided by R. Maddalena and 
T. Dame (priv. comm.). As with the other datasets, these radiation temperatures were 
gridded to the Mollweide projection with the same pixels as in the DIRBE maps and then 
smoothed by 3 pixels x 3 pixels. The negative tail of the histogram of the map of peak 
12 CO J = 1 — > radiation temperatures suggests a noise level of 0.15 K. This is consistent 
with the noise level in the map of integrated CO brightnesses for line widths of about 
2km-s _1 . As with the integrated CO map, background subtraction was not necessary. This 
map looks very similar to the integrated CO map and, consequently, is not shown. 
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3. Models and Results 

The ratio L,(240 /zm)/I( 13 CO), abbreviated as r 240 , is shown as a contour map in Figure 5 
and is plotted against the 140 //m/240 color temperature, or Td, in Figure 6. The points in 
Figure 6 represent the high signal-to-noise positions: those positions with signal at or above 
5-a in I„(140//m), I„(240//m), and I( 13 CO). The total number of points in the plot is 674, 
but these are not totally independent. Given that the effective beamsize is about 3 pixels x 
3 pixels, the effective number of points is about 70. (This is actually a conservative estimate 
to the effective number of points, because, even though the pixels within the effective beam 
are not completely independent, neither are they completely dependent.) This number is 
important for assigning a confidence level to the model fits. 

The physical models used to interpret Figure 6 include single- component LTE models, 
single-component non-LTE models, and two-component non-LTE models. Except where 
mentioned otherwise, the models are based on the following assumption: 

The only physical parameters that change from one line of sight to the next are the dust tem- 
perature, T d , and the gas kinetic temperature, T K , while maintaining a constant difference, 
AT = T d — T K . Other physical parameters such as gas density, dust-to-gas mass ratio, dust 
mass absorption coefficient, etcetera are assumed to be constant from position to position. 

This will be referred to as the basic assumption. Obviously, this assumption is not strictly 
correct, but is nevertheless roughly correct in describing the overall trends in the data. Thus 
the x 2 statistic that quantifies the goodness of the fit is not only measuring the scatter due 
to observational uncertainties, but also the scatter due to the real physical variations not 
accounted for in the basic assumption. The best fit of a certain class of models is found by 
adjusting a set of model parameters, aj, so as to minimize the x 2 ■ Because there are errors 
in both variables — i.e. x and y or, more specifically, T d and r 240 — the expression for x 2 
must reflect this: 
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and where (x oi ,y oi ) represents observed data point i, (x mi ({aj}),y mi ({aj}) is the point 
on the model curve closest to data point i, the {aj} are the parameters that characterize 
the model curve, and o-(x oi ) and a(y i) are the random errors for data point i in i and 
y, respectively. The closest point on the model curve to data point % is defined as that 
point on the curve that gives the minimum d\. Expressions (6) and (7) represent orthogonal 
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regression, in which the fit seeks to minimize the perpendicular distance between each point 
and the model curve. 

In the following sections, we examine the following cases: LTE, non-LTE one component, 
and non-LTE two components. 



In this section the 13 CO J = 1 — > line is assumed to be in LTE and optically thin. In 
practice, LTE means that the gas density is high enough that the line strength is not explicitly 
dependent on that density. The critical density of this line, at which the downward collisional 
transition rate is equal to the spontaneous radiative transition rate, is 3 x 10 3 cm~ 3 . Hence, 
this line is "close" to LTE for densities exceeding that critical density, or 10 4 -10 5 cm~ 3 or 
higher. Assuming LTE has the advantage that it is only necessary to fit one parameter, 
AT = T d - T K . The model curve to be fit is 



where Q is the rotational partition function for 13 CO and C BG is the correction for the cosmic 
background radiation. The reader is referred to Appendices A and B for details. Figure 7 
shows the model curves of equation (8) for AT = —16 to +16 K. 

The optimum value found for AT for the LTE models, i.e. AT = —4 ± IK, is far 
from satisfactory: as is seen in Figure 8, even though this curve goes through the center of 
the triangular cluster of points, it is systematically too high for the points with Td > 20 K. 
Indeed, the chi-square per degree of freedom is xl — 16.5. Fitting only to the points with 
T d > 20 K gives the result AT = +9 ± 1 K with xl = 9.0, but this curve is too low for the 
center of the triangular cluster of points. In addition, neither the AT = — 4 K curve nor the 
AT = +9 K curve match the high-r 240 or low-r 240 points in the triangular cluster. 

The uncertainties specified for the above AT value are the formal uncertainties due to 
the random errors. Normally, this uncertainty is determined by finding the values of AT 
for which the x 2 value is increased by 1. However, this is appropriate only if the chi-square 
per degree of freedom, xti * s near 1 ( see Press et al. 1992). This is not the case here: the 
xl is significantly greater than 1. Consequently, the uncertainty in AT was determined by 
finding the values of AT for which the x 2 was increased by the value xl- This is a more 
conservative estimate of the error in the fitted quantity. 



3.1. 



LTE Models and Results 
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Larger than the random uncertainties are the systematic uncertainties; these uncertain- 
ties include the calibration errors in the data and the uncertainties in the assumptions used 
in comparing data with theory. These uncertainties are listed below: 

1. The calibration of the observed r 240 or 1^(240 /im)/I( 13 CO) ratio. The calibration un- 
certainty of Ij,(240 (Km) is 12% (Hauser et al. 1998) and that of I( 13 CO) is less than 
20% (Nagahama et al. 1998). Combining the two uncertainties quadratically would 
imply a calibration uncertainty of 24% for r 240 . For simplicity, we adopt a calibration 
uncertainty of 20%. 

2. The uncertainty in the abundance of 13 CO, X( 13 CO). This affects the relationship 
between the CO column density and the H2 column density. Again, 20% is adopted as 
the contribution to the uncertainty in the model curve. 

3. The uncertainty in the dust optical depth to total gas column density, r„/N(H). This 
affects the relationship of 1^(240 /im) to N(H) (see equation B6). It is assumed that 
this contributes an uncertainty of 20% to the model curve. 

For simplicity, all of the above have adopted uncertainties of 20%, which is a number that 
is roughly appropriate for the above cases. Also, these uncertainties are independent of one 
another, so that they should be added quadratically. This gives an overall uncertainty of 
30-40% between the model curves and the data. 

Note that the systematic and random uncertainties are treated separately. The system- 
atic uncertainties represent a simple scaling up or down of the axes in plots of one quantity 
against another, while the overall distribution of points in a plot does not change. Hence, 
parameters determined from fits to the data, like line slopes, are simply scaled up or down 
accordingly. The random uncertainties, however, represent random reshapings of the distri- 
bution of points in a plot and must be included in the model fitting procedure. 

To test the effect of the systematic uncertainties on the result, the model fitting is 
repeated after applying scale factors to the model curves. The scale factors used in the new 
model fits deviate by as much as 40% from unity and are in steps of 20%: 0.6, 0.8, 1.0, 
1.2, and 1.4. As Figure 9 clearly shows, the total variation in the resultant AT values due 
to the systematic uncertainties is much larger than the error bars of the random or formal 
uncertainties. For the fits to all the high signal-to-noise positions, the AT varies from —24 K 
to +6 K. For the fits to the high signal-to-noise positions with Td > 20 K, the AT varies 
from —4 K to +9 K. Such a strong variation of AT, a range of 30 K for the fits to all the high 
signal-to-noise data, means that AT cannot be determined in this simple LTE case. This is 



-19- 



not surprising given that the fits are not very good. Also, since only one parameter varies 
anyway scaling the model curves varies that one parameter — AT in this case. 

Another test of the model fit is to compare the gas column density as derived from the 
continuum observations, N d (H), to that derived from the 13 CO J = 1 — > line, N 13 (HI + 
2H 2 ). The N d (H) values are determined from a trivial rearrangement of equation (B6) in 
Appendix B and the Ni 3 (HI + 2H 2 ) values come from 

N 13 (HI + 2H 2 ) = N(HI) + 2N 13 (H 2 ) , (9) 

in analogy to equation (B3). The Ni 3 (H 2 ) comes from equation (A27). The value of T K used 
in (A27) for each position is given by the value of T d for that position minus the AT of the 
model fit. As stated above, AT = — 4K, but AT = +9K was used for those positions with 
T d > 20 K. A plot of N 13 (HI + 2H 2 ) versus N d (H) is shown in Figure 10. The error bars 
represent the random uncertainties; the details of which are discussed in Appendices C and 
D. The systematic uncertainty in N d (H) is 30%. This includes the calibration uncertainty 
of 12% for 1^(240 jum), the calibration uncertainty of 6% for 1^(140 /xm)/I I/ (240/xm), and 
the uncertainty due to the background subtraction of 10%. The systematic uncertainty 
in N 13 (HI + 2H 2 ) is 20% due to the calibration uncertainty in I( 13 CO). The quadrature 
sum of these two uncertainties means that the systematic uncertainty in the slope in the 
Ni 3 (HI + 2H 2 ) versus N d (H) plot is about 40%. The slope of the plotted straight line and the 
y-intercept, as stated above, were determined from a routine that considers the uncertainties 
in both x and y, similar to the orthogonal regression method discussed previously. The 
uncertainties in the slope and y-intercept were determined by finding the changes in each 
that would raise the total x 2 by one, provided that xl — 1- Since xl was greater than 1, 
the uncertainties in the slope and y-intercept were determined after scaling up the errors in 
x and y by the same amount until xl — 1- The slope of the plot, 1.06±0.03, suggests that 
the column densities derived from the 13 CO line agree favorably with those determined from 
the continuum. However, the scatter is large: xl = 16.6. This is clearly a reflection of the 
bad fit of the single-component, LTE model to the data. 

As a comparison, Figure 10 also includes a plot of the gas column density as derived 
from applying the usual conversion factor to the integrated intensity of 12 CO J = 1 — > 0, 
I( 12 CO). As stated previously, the conversion factor employed was N(H 2 )/ J T H dv = 2.6 x 
10 20 cm~ 2 • (K ■ km ■ s -1 ) -1 . The slope in the fitted line is 1.40±0.04. This suggests that a 
more appropriate value for the conversion factor is 1.9 x 10 20 cm~ 2 ■ (K • km • s -1 ) -1 ; this will 
discussed further in Section 3.6. This value for the X-factor is only a crude average for the 
Orion clouds, because there is large scatter about the straight-line fit: xl — 28.3. Despite 
this large scatter, the 12 CO-derived column densities are within a factor of 2 of those derived 
from the continuum observations for about 80% of the points. So even on sub-cloud scales, 
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the X-factor seems to give column densities that are correct to within a factor of 2 for most 
positions. 

Comparison of the scatter in the upper and lower panels of Figure 10 may provide 
insights into the reliability of 12 CO-derived column densities. While the xl for the lower- 
panel straight-line fit with the 12 CO-derived column densities is nearly double that for the 
upper panel, the samples of points in the two panels are not the same. The sample with 
the 13 CO-derived column densities includes only those positions with greater than 5-cr in 
1^(140 (Km), I„(240//m), and I( 13 CO). This is 674 points. The sample with the 12 CO-derived 
column densities includes those positions with greater than 5-cr in I„(140 /im), I„(240 //m), 
and I( 12 CO). This is 1053 points. If the same 674-point sample is also used for the 12 CO- 
derived column densities, then xl = 17.5 for the straight-line fit in the bottom panel. This is 
comparable to the xl — 16.6 found for the straight-line fit to the upper panel with the 13 CO- 
derived column densities. At face value, the comparable xl values would suggest that the 
12 CO-derived column densities are almost as reliable as those derived from 13 CO. However, 
considering that the single-component LTE models used here do not properly characterize 
the data, this conclusion would be premature. 

Another important point to consider is that of the area filling factors of the gas. The 
model fitting has the potential of determining the gas kinetic temperature. Since the 
12 CO J = 1 — > line is optically thick, this line could also give the gas kinetic temper- 
ature, if we assume LTE. The peak radiation temperature of this line, T R ( 12 CO), would give 
T K . However, T R ( 12 CO)-derived T K values differ from the model T K values. This discrep- 
ancy is easily resolved by considering that the gas does not fill the beam within each narrow 
velocity interval. The radiation temperature would then be given by 

T» = r?» X(T K ) [1 - exp(-r,)] , (10) 

where 
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Expression (10) assumes that, at each velocity within the line profile, the telescope beam has 
an area devoid of gas emission and an area with gas emission. The t v is the optical depth 
of the gas averaged in the appropriate way over the gas-emitting area at frequency v within 
the line profile. The r) f (v) is the area filling factor and is the fraction of the beam area that 
has gas emission at the frequency v within the line profile. Given that the 12 CO J = 1 — > 
line is optically thick, equation (10) can be written as 
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Figure 11 shows the r] f (u) for the peak of the 12 CO J = 1 — > line profile determined from 
(12) plotted against the continuum- derived column densities. The gas kinetic temperatures 
were computed in the same way as for the column density calculations. There are two 
important things to notice in this plot. One is that, on average, positions with higher 
column densities have higher filling factors. This is expected: as the gas column density 
increases, there is an increasing number of cloud substructures or clumps along the line of 
sight. On average, more and more of the beam area is filled with clumps. Another thing 
to notice is that all the area filling factors are less than unity, as expected for a geometric 
filling factor. 

In short, the single-component, LTE models are not very successful in explaining the r 240 
versus T d plot. At best, these models only crudely reproduce the overall trend in the data 
and allow only a rough agreement between 13 CO-derived column densities and dust-derived 
column densities. 



3.2. Non-LTE Models and Results 

In this section, non-LTE effects on the 13 CO J = 1 — > line emission are considered. 
Also, the effects of non-negligible optical depths are also considered. To model the non- 
LTE line emission, the Large Velocity Gradient (LVG) method was used (Goldreich & Kwan 
1974). The advantage of this method is its simplicity. Photodissociation region (PDR) 
models (e.g., see Mochizuki & Nakagawa 2000; Storzer et al. 2000; Tielens & Hollenbach 
1985) are considerably more complicated and not necessary for the crude modeling done 
here. The LVG method assumes that every point in the gas is identical to every other point, 
with the exception of the velocities. In the spherical geometry case, every point in the gas 
is surrounded by points moving away from it in a spherically symmetric pattern with more 
distant points moving away even more rapidly — similar to a cosmological Hubble flow. 
Hence, every point in the gas is radiatively coupled to, and at the center of, a spherical 
region. Outside this region, the line emission of the gas is redshifted beyond that of the gas 
at the center of the sphere and, consequently, radiatively decoupled from it. This spherical 
geometry can then estimate the line emission from a spherical clump. This is only an 
estimate because the LVG method cannot accurately estimate the emission from the edges 
of a clump. Also, the velocity pattern within the clump would not be spherically symmetric. 
Nevertheless, the geometry of the region that is radiatively coupled to each point within the 
clump may be roughly spherical. Each clump is assumed to be radiatively decoupled from 
the other clumps; the other clumps that have radial velocities within the clump velocity 
width of a given clump are assumed to fill a small solid angle with respect to An steradians. 
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Once some basic parameters are chosen — the molecule, cosmic background tempera- 
ture, geometry, and the maximum number of levels to be included in the computation - 
three parameters are specified: the kinetic temperature, gas density, and a third parameter. 
This third parameter can be the abundance per unit velocity gradient, X( 13 CO)/(dv/dr) for 
13 CO, or the column density per velocity interval, N( 13 CO)/Av. It is common to choose 
X( 13 CO)/(dv/dr) as the third parameter, but this has a couple of problems. One problem is 
that, by itself, X( 13 CO)/(dv/dr) has little physical meaning. It says, in effect, that a region of 
the gas with a certain velocity width has a certain size (assuming fixed abundance); this says 
very little about the excitation by itself; the density must also be known. On the other hand, 
N( 13 CO)/Av states that a region of gas with a particular velocity width has a certain column 
density; this provides strong clues to the level of radiative trapping. For example, in the range 
of parameter space described in the next paragraph, given the optical depth of the J = 1 — > 
line of CO for given T K and N( 13 CO)/Av values, one can estimate the optical depth to within 
a factor of ten for other values of T K and N( 13 CO)/Av about 80% of the time even without 
knowing the gas density (provided that that unknown density is held constant). Another 
problem with using X( 13 CO)/(dv/dr) is that the effects of collisions and radiative trapping 
are mixed together, making interpretation of observational data more difficult. Specifically, 
the gas density, n(H 2 ), affects not only the collisions but also the column density. Because 
N( 13 CO)/Av = n(H 2 )X( 13 CO)/(dv/dr), raising n(H 2 ) for fixed X( 13 CO)/(dv/dr) also raises 
N( 13 CO)/Av, thereby raising the line optical depth while also raising the collision rate. This 
results in misleading impressions of the density sensitivity of line emission. For example, it 
has been mistakenly found that the 13 CO J = 1 — > 0/ 12 CO J = 1 — > intensity ratio is sensi- 
tive to n(H 2 ) (e.g., Sakamoto et al. 1997; Tosaki et al. 2002). This intensity ratio when plotted 
as a function of n(H 2 ) for fixed X( 13 CO)/(dv/dr) does show a strong variation with n(H 2 ), 
but this is because the rising n(H 2 ) also results in rising r( 13 CO J = 1 — > 0) and, of course, 
the 13 CO J = 1 -> 0/ 12 CO J = 1 -> intensity ratio is sensitive to r( 13 CO J = 1 -> 0). 
(In fact, r( 13 CO J = 1 — > 0) is equivalent to this intensity ratio in the LTE, optically thin 
limit.) Being unaware of this coupling between optical depth and density has led to errors 
in interpretation of the data. In one case, the authors went so far as to conclude that the 
molecular gas in the galaxy M 51 must have temperatures of several hundred kelvins on scales 
of hundreds of parsecs, based on a comparison of the J = 1 — > lines of HCN, 13 CO, and 
12 CO (Tosaki et al. 2002); this despite the much simpler and much more likely interpretation 
of low-temperature gas with optically thin 13 CO J = 1 — > emission. Therefore we choose 
N( 13 CO)/Av to be the third parameter used in the LVG models. 

The specific range of parameter space used in the modeling is as follows: 



• T K = 2.8 to 40.7 K in steps of 0.1K. 
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Zog>io[n(H 2 )(cm 3 )] = 1 to 5 in steps of 0.25. 



log 



10 



N( 13 CO) 
Av ' 



cm 



(km • s" 



15.5 to 18.3 in steps of 0.2. 



The reason for the chosen lower limit on N( 13 CO)/Av has to do with the observed col- 
umn densities of the Orion clouds. The majority of positions, 93%, have N(H I + 2H 2 ) < 
1 x 10 22 cm~ 2 . Negligible N(HI) and cloud line widths of 2 km • s _1 mean that this limit 
corresponds to N(H 2 )/Av = 2.5 x 10 21 cm~ 2 • (km • s -1 ) -1 . Using the abundances given in 
Appendix A, this is equivalent to N( 13 CO)/Av = 3.2 x 10 15 cm" 2 • (km • s -1 ) -1 . This column 
density per velocity interval represents a large-scale cloud property, for some positions in 
the Orion clouds, and must not be confused with the column density per velocity interval 
within the individual substructures or clumps. The N( 13 CO)/Av values within the clumps 
can be higher than the large-scale values for a cloud. However, the values for N( 13 CO)/Av 
cannot be much smaller than the large-scale value without also having many clumps within a 
given velocity interval that are radiatively coupled. Radiatively coupled clumps means that 
the effective N( 13 CO)/Av is roughly the sum of the N( 13 CO)/Av values for the individual 
clumps. Consequently, values much less than the large-scale cloud N( 13 CO)/Av value are 
not meaningful. Therefore, a lower limit of N( 13 CO)/Av = 3.2 x 10 15 cm~ 2 • (km • s" 1 ) -1 
is adopted, which is not too low for 93% of the high signal-to-noise positions (i.e. those 
positions represented in Fig 6) in the Orion clouds. 

To develop the expression for r 240 that uses LVG model results, we start with 

L = 2 B„(T d ) ^x d /i H m H N(H 2 ) , (13) 

where I„ is the continuum intensity. The above expression is a combination of equations (Bl), 
(B2), and (B3) with N(HI) = 0. The integrated strength of the 13 CO J = 1 -> line can be 
written as 

I( 13 CO) =T R ( 13 CO) Av , (14) 
which is exact if Av is defined appropriately. The quotient of (13) and (14) can be expressed 

N( 13 CO)" 



as 
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The subscript mod means that the quantity is either an LVG model input, [N( 13 CO)/A 
or output, [T R ( 13 CO)] mod . Using the values of the physical parameters from Appendices A 
and B, equation (15) becomes 
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Setting A to 240 fim gives 

r„ = 5.53 x 10- [exp(|?) - l] ["Ml [T^CO)]^ . (17) 

L mod 

Even though AT and n(H 2 ) do not explicitly appear in (17), the [T R ( 13 C0)] mO(i is an LVG 
model result, so it depends on T K (= T d - AT), N( 13 CO)/Av, and n(H 2 ). The effects of 
varying the parameters AT, N( 13 CO)/Av, and n(H 2 ) on the model curves are shown in 
Figures 12, 13, and 14. 

The best fitting curve is shown in Figure 15. The solid curve is fit to all data at more 
than 5-a in 1,(140 //m), 1^(240 /im), and I( 13 CO). The best fit parameters are AT = — 1±1 K, 
N( 13 CO)/Av= (3.2± 1 ^ 8 )xl0 15 cm- 2 -(km-s- 1 )- 1 ,andn(H 2 ) = (1.0± ? 4 ) x 10 5 cm" 3 . Note 
that it is not possible to put a firm lower limit on N( 13 CO)/Av nor a firm upper limit on 
n(H 2 ). In the former case, N( 13 CO)/Av is low enough that 13 CO J = 1 — > is close to the 
optically thin limit, and r 240 is insensitive to N( 13 CO)/Av in this limit. In the latter case, 
n(H 2 ) high enough that 13 CO J = 1 — > is close to LTE and r 240 is insensitive to n(H 2 ) in that 
limit. The chi-square per degree of freedom, xt * s 16.9, similar to that for the LTE model 
fit. The dashed curve is fit to the same 5-cr data, but only for Td > 20 K. The parameters 
for this curve are AT = -3 ± 4K, N( 13 CO)/Av = (3.2 ± \ 8 ) x 10 15 cm" 2 • (km • s" 1 )" 1 , and 
n(H 2 ) = (5.6 ± 25) x 10 3 cm -3 . The goodness of fit is xl = 9.98. Again, as for the LTE case, 
the LVG one-component model clearly does not well characterize the data. 

The poorness of the fit is also illustrated in Figure 16, which shows plots of the parameter 
variations with the scale factors. The total range in AT is 13 or 14 K. N( 13 CO)/Av varies 
by a factor of 7 for the solid curve and does not vary for the dashed curve. n(H 2 ) varies by 
nearly two orders of magnitude. If we ignore the cases with limits instead of error bars, then 
we see again that the systematic variations are larger than the random errors (except for 
the dashed curve showing the N( 13 CO)/Av variation). The particularly large uncertainty in 
n(H 2 ) is understandable because the 13 CO J = 1 — > line loses its density sensitivity for 
n(H 2 )>/eu; x 10 3 cm~ 3 . In short, as in the LTE case, it is difficult to determine the values 
of the parameters with certainty. 

Again the continuum-derived gas column density, Nd(H), is compared with that derived 
from the 13 CO J = 1 — > line, Ni3(HI + 2H 2 ). The Nd(H) values are determined as in 
the LTE case (see Section 3.1). The Ni 3 (HI + 2H 2 ) values come from equation (9), where 
Ni3(H 2 ) is derived from 

[N( 13 CO)/Av] d 13 
Nl3(H2) = XPCO) [T R p C0 J = l-0)] morf I( C0) ' (18) 

where the subscript mod refers to LVG model quantities (i.e. model input or output pa- 
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rameters), X( 13 CO) is the value adopted in Appendix A, and I( 13 CO) is the observed in- 
tensity. Each point on the model curve has particular values of [N( 13 CO)/Av] and 
[T R ( 13 CO J = 1 — > 0)] mod associated with it. Each data point was then associated with a 
point on the model curve. This was done in the same way that the x 2 values were computed: 
for each data point, the minimum of the distances to all the points on the model curve was 
the minimum of all values given by equation (7). (Also see discussion of this in Appendix E.) 
Note that the model curve used is the solid curve in Figure 15 for points with T d < 20 K and 
is the dashed curve in Figure 15 for points with T d > 20 K. The random errors associated 
with the Ni 3 (H 2 ) column densities were computed as described in Appendix E. The plot 
of N 13 (HI + 2H 2 ) versus N d (H) is shown in Figure 17. The straight-line fit to the points 
was done in same way as for Figure 10 in the LTE case. The slope of the line, 1.08 ± 0.02, 
suggests that there is good overall agreement between the continuum-derived and spectral 
line-derived column densities, except that the scatter is large: xt = 15.5. This is only 
marginally smaller than the xt = 16.6 obtained for the LTE case. Hence, the improvement 
from LTE, one-component models to LVG, one-component models is minor. 

The filling factors as a function of spectral line derived column densities in the LVG, one- 
component case are plotted in Figure 18. The filling factors were derived using equation (12). 
The T K used in that equation was determined from the observed T d at each position minus 
the appropriate AT from the model fits. (Even though these are non-LTE models, the 
12 CO J = 1 — > line is still essentially thermalized at these densities — i.e., n(H 2 ) = 1 x 10 5 
and 5.6 x 10 3 cm -3 — due to its high opacity. Consequently, the excitation temperature 
of this transition is close to T K and simply using T K in equation 12 is appropriate.) As in 
the LTE case, the filling factors, on average, increase with increasing column density. Again 
there is a saturation effect: the slope of the plot decreases with increasing column density. 

In short, the single-component, LVG models do not successfully explain the r 240 versus 
T d plot and are, at best, only marginally better than the single-component, LTE models. 



3.3. Two-Component, Non-LTE Models and Results 

As discussed in the previous sections, single-component models only poorly characterize 
the data. In this section we assume that there are two emission components. "Compo- 
nent 0" is present on all lines of sight and has constant properties on all those lines of sight. 
Specifically, the dust temperature, gas kinetic temperature, column density per velocity in- 
terval, and gas density are all constant. "Component 1" is similar except that the dust and 
gas temperatures can vary from line of sight to line of sight while maintaining a constant 
dust— gas temperature difference, AT, as stated in the basic assumption. The physical mo- 
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tivation for these components is straightforward. Component represents molecular gas at 
the surface of the Orion clouds and is heated only by the general interstellar radiation field 
(ISRF), giving it a constant temperature. Component 1 is the molecular gas deeper into 
the Orion clouds and is largely shielded from the general ISRF by component 0. Therefore, 
those positions in component 1 without stars will have colder dust temperatures than in 
component 0. In contrast, those positions in component 1 with stars will have warmer dust 
than in component 0. In this simple picture, component has a constant temperature and 
component 1 has a varying temperature from position to position, varying above and below 
the dust temperature of component 0. 

We now need expressions relating the observed emission to the physical parameters 
of the gas and dust in order to construct the two-component models. The expressions 
derived here will be necessary for the simulations described in Paper II. Construction of the 
models requires certain simplifications. For example, it is assumed that component has 
the same AT as component 1. Also, as discussed in Section 3.2, the clumps are assumed to 
be radiatively de-coupled from one another. Specifically, for these two-component models, 
it is assumed that the clumps of component are de-coupled from those of component 1. 
Therefore the emission of the two components can simply be added together. Consequently, 
the specific intensity of the continuum emission averaged over the beam, I„, is given by 

L = giN cl f„(T dl ) + g N c0 f,(T d o) , (19) 

where f„(T d ) is defined as 

( v\P 

f„(T d ) = K u0 x d fi H m H — B„(T d ) , (20) 

combining equations (Bl), (B2), (B5) and removing the column density factor. N ci is the 
column density of hydrogen nuclei (i.e. 2N(H 2 ) for purely molecular hydrogen) of each clump 
of component 1. gi is the number of clumps of component 1 on each line of sight averaged 
over the beam area. N c o and go are the corresponding quantities for component 0. The total 
number of clumps per line of sight averaged over the beam from both components, g, is given 
by 

(21) 



aiN c i + a N c0 

where N is the observed column density (i.e. the column density averaged over the beam) and 

ai, a are the fractions of clumps belonging to components 1 and 0, respectively. Obviously, 
then 

a + ai = 1 . (22) 

It follows then that 

gl = aig= N cl + c N c0 (23) 
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and 



where 



go = a g 



c N 
N c i + c N 



cO 



a 

c = — 

ai 



(24) 



(25) 



Notice that g = gi + g follows naturally from (21), (22), (23), and (24). Equation (19) can 
now be written as 



I, 



N 



N c i + c N c0 



N cl f„(T dl ) + 



N 



Nci + c N c0 



c N c0 f I/ (T d0 ) 



(26) 



To compare the continuum emission strength with that of a spectral line, the cloud's velocity 
structure must be considered. Even though the dust grain velocities have no appreciable 
effect on the continuum emission, we can imagine examining the emission of only those dust 
grains associated with the 13 CO molecules emitting in a certain part of the 13 CO J = 1 — > 
line. As such, the right side of (26) is multiplied by where Av is the observed line 

width and Av c is the clump line width. For simplicity it is assumed that all clumps have 
the same velocity width. Consequently, (26) can be written as 



I* = Ci 
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Av r 
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(27) 



(28) 



Analogous to gi, Ci is the number of clumps of component 1 per line of sight averaged over 
the beam within a velocity interval equal to a clump velocity width. The parameter Co (which 
is not analogous to g ) is the ratio of the number of clumps of component to those of 
component 1. (Note that clump velocity widths of the two components are assumed to be 
identical. Generalizing to the case of different clump velocity widths of the two components 
is trivial.) Analogous to the continuum emission equation (i.e. equation 27), the equation 
for the integrated radiation temperature of the 13 CO J = 1 — > line is 



I( 13 CO) = ci 



T R1 (T, 



Nci 
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Av 
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(29) 



where T R1 and T R0 are the component-1 and component-0 radiation temperatures, respec- 
;i and T K0 are the gas kinetic temperatures in components 1 and 0, respectively, and 



tively, T 

n c i and n c o are the molecular gas densities in the component-1 and component-0 clumps, 
respectively. Note that |f is and gf is The LVG results, T R1 and T R0 , 
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depend directly on Ncl ^ v CQ - ) and Nc0 ^ v CQ - ) , respectively. Consequently, T R1 and T R0 as writ- 
ten in equation (29) also implicitly depend on the 13 CO to H 2 abundance ratio, X( 13 CO). 
Dividing (27) by (29) gives the desired r 240 : 

fef,(T dl ) + c ^(T d0 ) 

(30) 



240 T R1 (T dl - AT, £^,n cl ) + c T R0 (T d0 - AT, |^,n c0 ) 
where the following equations were used: 

T K1 = T d i — AT (31) 

and 

T K0 = T d0 — AT . (32) 

Again, the same AT is assumed for both components. 

Equation (30) is used to determine the r 240 ratio for the two-component, non-LTE case. 
The T R1 and T H0 values are determined from the LVG method. These models have seven 
input parameters: AT, Co, T d0 , n c0 , ■^ L , and n cl . The component-1 dust temperature, 
T dl , varies, resulting in model curves of r 240 versus dust temperature. The dust temperature 
used in the plots of the model curves is the 140 //m/240 yum color temperature, T dc . To solve 
for T dc , the expression (27) is evaluated at vg corresponding to Ag = 140 /im and at v\q 
corresponding to Ai = 240 /im and these are divided to yield: 

= feUT dl )K^(T dl ,/?) + co^f^CTdo^CTdo,/?) 
^ fi/io(Tdi)K J/10 (T dl , P) + c ^f, 10 (T d0 )K, 10 (T d0 ,/5) 

which uses the same notation as expression (4) in Section 2.1. The color corrections, 
Ki,(T d ,/3), are used in order to correctly assess the relative contributions of the two compo- 
nents to DIRBE's bands 9 and 10. Equation (4), with T dc in place of T d , is then solved for 
T dc . This allows a valid comparison of the model T dc values with the DIRBE observed T dc 
values. 

The master grid of T K and n(H 2 ) values is the same as that for the one-component, 
LVG models. The master grid values of log w N( - A ^°- > ( cm" 2 ■ (km ■ s" 1 ) -1 ) is extended 
down from 15.5 to 14.3. This does not necessarily violate the constraint that the minimum 
column density per velocity interval be about that of the cloud on the scale of a beam, 
because we are now dealing with two components. Having the combined column density per 
velocity interval of the two components above this minimum is still required. 

The best-fit curve is shown in Figure 19 and is fitted to all the data at more than 5-<r in 
1^(140 pLm), 1^(240 (Km), and I( 13 CO), and the parameters are listed in Table 1. The curve is 
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characterized by the triangular loop that exists between T dc ~ 15 and ~ 19 K and r 240 ~ 7 
and ~ 40 MJy-sr~ 1 -(K-km-s -1 ) -1 . This loop is a consequence of the two components emitting 
in the Wien limit. The formal uncertainties are not listed because they are extremely small 
and not really meaningful when compared with the systematic uncertainties. The formal 
uncertainties for AT and Tdo are less than 3 x 10~ 4 K and for the other parameters the 
formal relative uncertainties are less than 3 x 10~ 4 . Unlike the one-component model curve, 
the two-component model curve is a reasonably good fit to the data, with xl — 5.7. The 
F-test states that the two-component model is better than the one-component models at the 
99.9% confidence level. (Note that the number of degrees of freedom adopted for estimating 
this confidence level was about 1/9 the number in Table 1, because the points in Figure 19, 
which represent individual map pixels, are not independent. Given that the number of 
pixels per effective beam is about 9, the effective number of degrees of freedom is reduced 
accordingly. This reduction results in a lower confidence level than would have occurred 
without this reduction.) Figure 19 explains the overall trend in the data and even explains 
a large part of the triangular cluster of points (see Section 3). Nevertheless, this fit still has 
certain problems: the curve does not follow the curvature in the data for T dc > 20 K and it 
does not account for the low-r 240 points in the triangular cluster where T dc > 18 K. While the 
best-fit curve does not explain these low-r 240 points, some combination of parameter values 
could account for these points: if the component-0 parameters were allowed to vary from 
position to position (i.e., ignoring the basic assumption), for example. Also, allowing Tdo to 
vary between about 18 K and 20 K could fit these points. Nevertheless, given the simplicity 
of these two-component models, the current fit is reasonable. In addition, as mentioned in 
Section 3.1, these points could also be explained by uncertainties in the observations. 

Even though the curve in Figure 19 is an adequate fit, the parameters in Table 1 still have 
appreciable uncertainty. Figure 20 shows the systematic effects on the resultant parameter 
values when the scale factor is changed. Parameters like AT and T^o are known extremely 
well, ±1 K or better, while parameters like Nc0 ^ v CQ - > and n c o are uncertain by 2 or 3 orders 

of magnitude. Still others like c , Ncl ^ v CO ' ) , and n c \ are apparently known to within about 
an order of magnitude or better. In the case of the one-component models, ascertaining any 
of the parameter values was difficult. In the case of the two-component models, again some 
parameter values are very uncertain, but some are apparently very well known indeed. 

The systematic uncertainties in the resultant parameter values were also tested by vary- 
ing the starting search grid. There are enough parameters in the two-component models that 
a search grid with both sufficient range in all parameters and sufficient resolution (i.e. grid 
spacing) in all parameters was not possible. Therefore, the starting search grid had a large 
range in all parameters and poor resolution (i.e. large grid spacing). After the initial search, 
a new grid with better resolution and smaller range moved closer to the final solution. This 
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was repeated a few times to home in on the final solution. In total there were 3 to 5 search 
grids (including the starting grid) for each final solution. Of course, one wonders about the 
final solution's stability: does the starting search grid affect the final answer? In the ideal 
case the answer is no. Nevertheless, this was tested by adjusting the starting search grid's 
position to see the effect on the final answer. In total, four different starting search grids 
were tested and the results were examined. The results did indeed change depending on the 
starting grid. Nevertheless, the ranges in the results' variations had nearly identical sizes (in 
a linear sense for AT and T d0 and in a logarithmic sense for the other parameters) to those 
found from varying the scale factors. 

As in the previous sections, the gas column density as derived from the continuum 
observations, N d (H), is compared with that derived from the 13 CO J = 1 — > line, N 13 (HI + 
2H 2 ). The Ni 3 (H 2 ) values are derived from a modified version of equation (18) that was used 
for the one-component, LVG model: 



N c i( 13 CO) 
Av c 



+ C 



Ncu( 13 CO) 
Av c 



N 13 (H 2 ) = - r- I( 13 CO) . 

XPCO) {[T R1 PC0 J = 1 - 0)] mod + c [T R0 (-CO J = 1 - 0)] mod } 

(34) 

Using this equation is similar to using equation (18) for the one-component, LVG models. 
As for those one-component models, the random errors associated with the N 13 (H 2 ) column 
densities were computed as described in Appendix E. The N d (H) values must now be com- 
puted and require corrections from the one-component cases of the previous sections. This 
correction is 

Nd(H) N,m , . 

Ndi(H) N 1M ' [6b) 

where N d (H) is the "true" dust-derived column density inferred from the observations in 
the two-component case, N dl (H) is the dust-derived column density that would be inferred 
from the observations assuming one component only, N 2 m is the model column density 
for a two-component model, and Nim is the model column density that would result if a 
single component were assumed. Specifically, this single component would have the color 
temperature, T dc , and the specific intensity of the two-component model, as given by (27). 
So the right side of (35) is the correction factor that would be applied to the model column 
density to obtain the "true" model column density. Equation (35) then applies this correction 
factor to the actual observations. If all the column densities in (35) are beam-averaged 
column densities, then 

N 2M = gi(N cl + c N c0 ) (36) 

and 

Nim = giN c . (37) 
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Equation (36) is combination of equations (21), (23), (24), and (25). The N c in equation (37) 
is the effective column density of a single clump for the one-component "version" of the two- 
component model. It is given by the expression, 



N c 

I^o = Ci^f i , 10 (T dc )Av 



(38) 



which is the one-component analog of expression (27) that uses u w in place of v. Equating 
(38) with (27) and rearranging yields 
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Expression (39) is substituted into (37), and (37) and (36) are substituted into (35). After 
this, a factor of X( 13 CO) is removed from both numerator and denominator to yield 
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(40) 



As done in the previous section for the one-component, LVG models, each data point was 
associated with the point on the model curve closest to it. The parameter values associated 
with that point on the model curve were then used in (34) and (40) to give the column 
densities N 13 (H 2 ), and therefore N 13 (HI + 2H 2 ), and N d (H) corresponding to that data 
point. The uncertainties in N d (H) are discussed in Appendix F. 

Figure 21 is the plot of the 13 CO-derived column densities versus the dust-derived gas 
column densities. Clearly the agreement is much better than for the one-component models: 
xl = 7.4. The slope of the straight-line fit, 1.08±0.01, is within the error of the corresponding 
slope for the one-component, LVG models. It is clear that the two- component models are 
significantly better than the one-component models. 

The filling factors as a function of spectral line derived column densities in the LVG, 
two-component case are plotted in Figure 22. The filling factors were derived using a 
two-component analog of equation (12). Specifically, an effective excitation temperature, 
Jv{eff)-> is determined from 



Meff) = 



X(T xl ) + c X(T x0 ) 
" 1 + c 



(41) 



where T xl and T x0 are the excitation temperatures of the 12 CO J = 1 — > line in com- 
ponents one and zero, respectively. This effective excitation temperature is the brightness 
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temperature that would be observed at frequency v if the gas filled the beam at the veloc- 
ity in the line profile corresponding to frequency v. (Note that this assumes that the line 
is optically thick. Note also that J v [ej f) includes the Rayleigh- Jeans correction and the 
correction for the background.) The excitation temperatures, T xl and T x0 , were computed 
with the following method. Each observed position has an associated point on the model 
curve. The point on the model curve has associated Tdi, Tdo, and AT values. These allow 
computation of the T K1 and T K0 values. The point on the model curve is also associated 
with the other model parameters, which permit determination of the T xl and T x0 values 
from LVG models of the 12 CO line strengths. In particular, the J = 1 — > line of 12 CO is 
optically thick and, as such, radiative trapping keeps the line closer to LTE. Therefore, the 
conversion from kinetic temperature to excitation temperature is close to unity. Specifically, 
the following approximations are accurate to within a few percent: 

Txo - T K0 (42a) 

and 



T xl ~ T 
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(42b) 



Over the observed temperature range, the conversion of T K1 to T xl deviates from unity by 
no more than about an 8%. The filling factor is then given by 

"•<^m ' (43) 

The value of T) f (v) at the line peak is plotted in Figure 22. As in the one-component cases, the 
filling factors, on average, increase with increasing column density. Again there is saturation 
with increasing column density. However, for these models, the area-filling factors can have 
non-physical values (i.e. > 1). In fact, some observed positions seem to have T] f (v) up to 
about 6. Since these are area-filling factors, they cannot be greater than unity. Even if we 
were to reinterpret T] f (v) as the number of clumps per line of sight within a clump velocity 
width, there would still be a problem: the effective column density per velocity interval 
would be higher than, and therefore inconsistent with, the values of Ncl ( CQ ) anc l Nc °( _co) 
determined from the model fits. Only for r] f (v) < 1 is there no such problem, because some 
lines of sight within the beam, and within a clump velocity interval, can be free of clumps 
while others can have a single clump (with the properties given by the model fit). When 
rj f {u) > 1, lines of sight that have clumps will, on average, have more than one clump in 
a clump velocity interval, thereby forcing the effective N( 13 CO)/Av to be higher than the 
fitted Ncl ^ CQ - > and Nc0 ^ CQ - > values. In addition, given that this is a filling factor within a 
clump velocity interval and that 12 CO J = 1 — > is optically thick, detecting more than one 



-33- 



clump per line of sight within a narrow velocity interval is not possible, because the clump 
in front would block the emission of the ones behind. 

Nevertheless, the problem of the unphysical r] f (u) values has a solution. Both the nu- 
merator and denominator of (30) depend, explicitly or implicitly, on the product c |^-, and 
not on the individual Co and values alone, over a large range parameter space. Accord- 
ingly, over that range of parameter space, r 240 depends on that product and determining the 
parameters Cq and separately is very difficult. Indeed, Figure 20 shows that Cq and ^ L 
are roughly anticorrelated. The plot of their product shows that it is constant to within a 
factor of roughly 3 about some geometric mean. Consequently, restricting Co to some desired 
range can still yield physically reasonable fits. This will be examined in the next section. 

In summary, the two-component, LVG models are reasonably good in explaining the r 240 
versus T dc plot, which is in stark contrast to the one-component models. Examination of the 
systematic uncertainties suggests that AT and T d0 are known to within 1 K or better and 
n c i is known to within a factor of 2 about some geometric mean. Other parameters like Co 
and are known to within factors of 4 or 5 about a geometric mean, while the remaining 
component-0 parameters, and n c0 can vary by factors of ~ 30 about the geometric 
mean. While a different class of two-component models could possibly explain the entire 
triangular cluster of points in the r 240 versus T dc plot, the two-component models described 
here only explain the triangular cluster for T dc < 18 K. The points in the triangular cluster 
with T dc > 19 K are not accounted for with the current models. Another problem with the 
current models is that they predict area-filling factors that are unphysical for many of the 
observed positions. And this problem will be addressed in the following section. 



3.4. Two-Component, Two-Subsample, Non-LTE Models and Results 

As mentioned in the previous section, the two-component models suffer from a non- 
trivial physical drawback: the area-filling factors for some positions are greater than unity. 
Another problem is that the two-component models do not fit the T dc > 20 K points well. 
As discussed in that section, keeping c large enough will ensure that the filling factors 
have more reasonable values. Therefore, in this section we will consider models where Co 
is restricted to be > 1. To better fit T dc > 20 K points, we divide the points into two 
subsamples: one subsample with T dc < 20 K and the other with T dc > 20 K. We fit the two- 
component, LVG models with c > 1 only to the T dc < 20 K subsample. The T dc > 20 K 
subsample has a relatively simply trend and we need only fit the one-component models 
that were used in Section 3.2. Notice that we are relaxing the basic assumption somewhat: 
we are allowing parameters, in addition to T d and T K , to change from one subsample to 
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the other subsample. However, within each subsample we still impose the basic assumption. 
The parameter values resulting from the fits to the two subsamples are listed in Table 2. 
The resultant model curves are shown in Figure 23. 

Casual inspection of the curves in Figure 23 suggests that these models better charac- 
terize the data than the simple two-component models of the previous section. However, the 
appropriate combination of the xl values in Table 2 into a reduced chi-square that charac- 
terizes both curves together yields xl = 5.73. Table 1 shows that the two-component models 
of the previous section yielded a best fit of xl = 5.69. Hence, the model curves for the 
two subsamples fit the data with essentially the same goodness of fit as for the simple two- 
component models. To be consistent with the two-component model fit to the T dc < 20 K 
subsample, we could try to fit a two-component model to the Td c > 20 K subsample. This 
two-component fit to the T dc < 20 K subsample, when combined with that for the T dc > 20 K 
subsample, gives an overall xl — 5.34. Again, this is not significantly better than the sim- 
ple two-component models for the entire sample of points. Also, this model requires fixing 
T d0 at 18 K to get a reasonable fit. Therefore, only the one-component model fit to the 
Tdc > 20 K subsample is listed in Table 2 and is used in the column density and filling factor 
determinations. 

If we consider only that the goodness of fit is no better and that extra model parameters 
are required, then the models discussed in this section should be ignored in favor of the 
two-component models for the entire sample of points. However, the two-component, two- 
subsample model described here has three important advantages. One is that the model fits 
the T dc > 20 K points better. The second advantage is that the agreement between the dust- 
derived column densities, N d (H), and the 13 CO-derived column densities, N 13 (HI + 2H 2 ), is 
noticeably improved. Figure 24 shows the plot of Ni 3 (HI + 2H 2 ) versus N d (H). The xl of 
the straight line fit to the points in this plot is 5.3, compared with 7.4 in the corresponding 
plot for the simple two-component models. This an improvement with a 90% confidence 
level, according to the F-test. 

Another advantage to the two-component, two-subsample models is seen in Figure 25, 
which shows the area-filling factors, r] f (u), versus the 13 CO-derived molecular gas column 
densities, Ni 3 (H 2 ). The filling factors are now all (or almost all) less than or equal to unity, 
as desired. As was seen for the one-component models, the filling factors here again show 
saturation with increasing column density. 

The uncertainties are estimated from the work of the previous sections. For the T dc < 
20 K subsample, the uncertainties of the two-component model parameters are largely sys- 
tematic. As such, the uncertainties should be similar to those for the two-component models 
fitted to the entire sample of points (see previous section). However, our estimates of those 
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uncertainties for the component-1 density, n c \ must be modified. The previous section sug- 
gests that n c i is known to within a factor of 2 about a geometric mean of possible densities. 
However, the n c \ for the current section is 2 orders of magnitude larger than for the previous 
section. Hence, we can only place a lower limit of ~ few x 10 2 cm~ 3 on n c \. We obtain a 
higher density value in the current section, possibly because the densities in the T dc < 20 K 
subsample are 2 orders of magnitude higher than some kind of average density for the whole 
sample. Indeed, we find that the T dc > 20 K subsample has a density that is lower by more 
than an order of magnitude. However, the T dc < 20 K subsample has always dominated the 
results of the fits; the results obtained for this subsample are nearly identical to those for the 
whole sample. Therefore, this method likely places only a rough lower limit on the density. 
This is expected because the 13 CO J = 1 — > transition is insensitive to densities higher 
than about 10 3 cm -3 . Thus the uncertainties in the parameter values for the T dc < 20 K 
subsample are the same as those for the entire sample points (see the previous section) with 
the exception that the component-1 density only has a lower limit of ~ few x 10 2 cm -3 . 

The parameter uncertainties for the Td c > 20 K subsample are again dominated by 
the systematics. The AT versus scale factor for the one-component models, i.e. Figure 16, 
shows that the fit to the T dc > 20 K points varies from —10 to +4K. Hence we can say 
that AT = —3 ± 7K. This uncertainty is only slightly larger (i.e. less than a factor of 2) 
than the formal uncertainty given in Table 2. The plot of N( 13 CO)/Av versus scale factor 
suggests that only the formal uncertainties in N( 13 CO)/Av are important. However, the two- 
component model fitted to the T dc > 20 K sample suggests larger systematic uncertainties. 
The Ncl ^ CO - ) value found (which corresponds to the N( 13 CO)/Av value in the one-component 
model) is less than a factor of 2 greater than the N( 13 CO)/Av value found for the one- 
component model. The systematic uncertainties (see the plot of Ncl ^ v CQ - > versus scale factor 
in Figure 20) suggest that N( 13 CO)/Av varies between about 10 15 and 10 16 cm~ 2 -(km-s -1 ) -1 . 
(Recall that it cannot be lower than about 10 15 cm -2 • (km • s -1 ) -1 because the large-scale 
observed N( 13 CO)/Av represents a lower limit.) As for the systematic uncertainty of the 
density, Figure 16 suggests that, again, only a lower limit is applicable. This lower limit 
is ~ 10 3 cm~ 3 . It is interesting that the two-component model results for the T dc > 20 K 
subsample suggest a component-1 density identical to that of the one-component model fit. 

One noticeable problem with the current models is the abrupt transition in the solid 
curve when crossing the T dc = 20 K boundary. A more natural-looking curve would have a 
smooth transition across this boundary. Examination of the different parameter values for 
both subsamples suggests that gradually decreasing the ^ value by factors of 4 to 6 when 
going from the T dc < 20 K subsample to the T dc > 20 K subsample (thereby becoming the 
N( 13 CO)/Av value for the one-component model for the latter subsample) could possibly 
produce the desired smooth curve. Note that other solutions, such as varying n c i would not 
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produce a satisfactory solution, because it would not merge well with the low-Td c curve. Also, 
the component-0 parameters Cq, and n c o have little affect on the overall curve shape 
and are irrelevant for the Td c > 20 K subsample. Therefore, the Td c > 20 K subsample 
clumps are likely to have N( 13 CO)/Av values that are factors of a few lower than those in 
the Tdc < 20 K subsample. 

The dashed curve on the Td c > 20 K side of Figure 23 suggests another possibility. The 
current models do not account for the vertical spread of the points (neither in the Td c > 20 K 
side nor in the triangular cluster mostly on the Td c < 20 K side) in the r 240 versus Td c plot. 
Also, they do not account for the horizontal spread of the triangular cluster, although this 
spread is probably due to the large uncertainties in Td c for some of these points. The 
dashed curve on the T dc > 20 K side passes through the high-r 240 . Accordingly, two or more 
subsamples could indeed account for these spreads if the curves for all the subsamples were 
applied to the whole range of observed color temperatures. However, what is the optimum 
way of separating these two (or more) subsamples? Separating them according to dust 
temperature is simplistic: the Td c > 20 K solution, for example, does not pass through points 
on the Tdc < 20 K side (i.e. the dashed line on that side) that are far from the Td c < 20 K 
solution (the solid curve on that side). The optimum separation of these subsamples would 
involve two (or more) curves for the whole color temperature range of that would correctly 
characterize all the points in the plot. In other words, the basic assumption in its simplest 
form is not correct. Nonetheless, the basic assumption imposed within each subsample would 
simplify the modeling. For now, we keep the current two-component, two-subsample model 
parameters because they represent an adequate characterization of the data. 

In summary, the two-component, two-subsample models produce a similarly good fit 
to the r 240 versus T dc plot as the simple two-component models, but fix the problems of 
the unphysical filling factors and the poor fit to the high-T d points. Also, the dust-derived 
and 13 CO-derived column densities agree better for these two-subsample models than for the 
models fitted to the whole sample. Also, the current two-component, two-subsample models 
hint that the optimum choice of subsamples has yet to be found. This will be left to future 
work. 



3.5. The Fiducial Points 

As a rough illustration of the limitations of the basic assumption and of the modeling, 
a few reference or fiducial points were chosen in the r 240 versus Td c plots. These points 
roughly follow the overall shape of the sample in these plots. Specifically, four of the fiducial 
points approximately represent the corners and centroid of the triangular cluster of points 
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and the other three fiducial points represent the overall trend in the high-T dc branch. The 
fiducial points are represented as diamonds in the r 240 versus Td c plots (i.e., Figures 6, 7, 8, 
12, 13, 14, 15, 19, and 23). Table 3 lists the fiducial points (in order of increasing Td c ) and 
the physical properties that pertain to each. Because of the basic assumption, the derived 
physical properties would be the same for all the points (except T K and Td). However, when 
a high-Tdc solution was found, as for all but the LVG two-component case, the physical 
conditions derived for the four fiducial points with T dc > 20 K were those of the high-T dc 
solution. The table lists these physical conditions for the four cases of the previous sections: 

Case 1: LTE, one-component models. 

Case 2: LVG, one-component models. 

Case 3: LVG, two-component models. 

Case 4: LVG, two-component, two-subsample models. 

Despite representing overall trends, the fiducial points do not necessarily represent the 
sample as a whole: e.g., the second, third, and fifth points as listed in Table 3, which have 
distances of 3 or more sigmas from the best-fitting model curves. The deviations of these 
three points from each model curve are not just statistical uncertainties; the models used 
here simply cannot account for the total vertical spread of the points in the sample. The 
basic assumption does not allow models that would "fill" the space occupied by the sample of 
points in the r 240 versus T dc plot. One alternative is to have two or three of these two-compt 
model curves fit to all the sample simultaneously. In practice, this is time-consuming and is 
left to a later date. 



3.6. Masses 

The derived column densities now allow us to estimate the masses of the Orion clouds. 
Maps of the column densities are presented in Figures 26 to 28. Figure 26 shows the column 
density maps derived from the observations using the LVG, one-component model results, 
Figure 27 the analogous maps for the two-component, two-subsample case, and Figure 28 the 
maps of column density for just component 1 on those lines of sight where this component has 
a temperature of 10 K or less (i.e., Tdi = T K1 < 10 K). The maps of column densities for the 
LTE case are not shown, because they closely resemble the maps of the LVG, one-component 
case. The upper panel in each of these figures shows the continuum-derived (dust-derived) 
gas column densities and each lower panel shows the 13 CO-derived gas column densities. 
Summing these column densities appropriately and using the adopted distance of 450 pc gives 
the masses (not including Helium) listed in Tables 4, 5, 6, and 7. The random uncertainties 
in these masses are small: relative uncertainties of the order of 0.5% for the continuum- 
derived masses, 0.25% for the 12 CO-derived masses, and 6.5% for the 13 CO-derived masses 
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(although the smaller Orion Nebula and NGC 2024 Fields have relative uncertainties of 20% 
for the 13 CO-derived masses). The systematic uncertainties (i.e., the uncertainties due to 
calibration and background-subtraction uncertainties) are more appreciable: 30% for the 
continuum-derived masses and 20% for the 13 CO-derived masses. The model results in each 
of the four cases were combined with the 13 CO J = 1 — > and H I observations to yield the 
column densities, as described previously, and these were summed to yield the masses in the 
tables. 

All the 13 CO-derived column density maps show a plateau with steep edges. This plateau 
is roughly 20° wide in longitude and about 6° in latitude. For example, one of the edges of 
this plateau can be seen running from galactic longitude 216° to 208° at latitude —10.25°. 
These are approximately the mapped area in 13 CO. Specifically, these edges occur just inside 
the mapped area's edges, where the observed 13 CO J = 1 — > line strength falls to zero. 
These plateau edges are also partly visible in the upper panel of Figure 27, even though this 
panel does not have the 13 CO-derived column densities. This is because these "observed" 
column densities depend on the two-component, two-subsample model results, which, in 
turn, require knowledge of the 13 CO line strengths. This model dependence is especially 
true of the cold gas column density maps in Figure 27, resulting in the clear visibility of 
the plateau edges in the northward and eastward (i.e. towards the top and left) directions. 
The southward and westward edges are the arbitrary boundaries imposed by selecting only 
those lines of sight with T d i < 10 K; column densities of component 1 on lines of sight with 
T dl > 10 K are set to zero. 

There are apparent discrepancies between the derived masses (i.e. Tables 4 and 5) 
and the column density versus column density plots (i.e., Figures 10, 17, 21, and 24). For 
example, the ratio of the 12 CO-derived masses in Table 4 compared to the continuum- derived 
masses of Cases 1 & 2 in Table 5 is 1.08, whereas the linear fit in the lower panel of Figure 10 
suggests a ratio of 1.40. There are two reasons for the discrepancy that are related to the 
signal-to-noise ratio of the observations. One reason is that the plot of column densities 
in Figure 10 only considers positions with high signal-to-noise data (i.e. more than 5-a in 
12 CO J = 1 — > and in the continuum bands). If all the positions in the Orion fields were 
included in the plot, the slope would change to about 1.3. The rest of the discrepancy is 
because of the difference between the slope of a line from linear regression and the ratio 
of the averages of the ordinate values and abscissa values. The former represents a sum of 
values that are weighted by 1/a 2 and divided by another sum of similarly weighted values. 
Determining the total masses by simply summing the column densities in the fields and 
dividing by another similarly determined sum represents the latter. In other words, the 
masses in the tables are not determined from weighted sums, whereas the slopes in the 
column density plots are indeed determined from such sums. Therefore, the slopes in the 
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column density plots are biased toward the points with higher signal-to-noise ratios even 
when these plots include the low signal-to-noise ratio points. This discrepancy in the ratio 
between the tables and the plots is simply a consequence of the scatter of the points in 
these plots. If there were less scatter, then the tables and the plots would give similar 
ratios for the 12 CO-derived over continuum-derived masses (assuming similar coverage in 
both types of maps). This affects the accuracy of our X-factor estimates. The slope from 
the lower panel of Figure 10 suggests an X-factor of 1.9 x 10 20 cur 2 • (K • km ■ s -1 ) -1 , whereas 
comparing the total masses in the tables for the one-component cases suggests an X-factor 
of 2.4 x 10 20 cm -2 • (K • km ■ s -1 ) -1 . Combining these two numbers gives an estimate of the 
X-factor of (2.1 ± 0.3) x 10 20 cm" 2 • (K • km • s -1 ) -1 . This is an uncertainty of about ±15%. 
Thus specifying the X-factor to greater accuracy (i.e., a smaller uncertainty) is meaningless. 

An interesting, but crude, comparison can be made with the Orion X-factor estimated 
from Figure 8 of Dame et al. (2001). In that figure, the predicted I(CO) profile that cuts 
along the Orion region is on average about 20-30% higher than that of the observed profile. 
Accordingly, the predicted X-factor of 1.8 x 10 20 cm~ 2 • (K • km • s -1 ) -1 for the entire Galaxy 
must be increased by factors of about 1.2 to 1.3 to give the average X-factor for Orion. This 
gives an average X-factor for Orion of about 2.2-2.3 x 10 20 cm~ 2 • (K • km • s -1 ) -1 . This is 
consistent with the value derived above. 

There is another apparent discrepancy between the plots and the tables; the total masses 
derived from 13 CO (i.e., in Table 4) disagree with those derived from the continuum (i.e., in 
Table 5). In this case, the ratios derived from the tables are more extreme (i.e., further from 
1) than from the column density plots. The problem here is that the 13 CO map does not 
cover the full area covered by the continuum maps. This results in two problems: a problem 
with the continuum-derived gas masses and a problem with the 13 CO-derived gas masses. 
The former problem is that the two-component-model results cannot give the continuum- 
derived gas column densities outside the area mapped in 13 CO. Specifically, an observed 
data point in the r 240 versus Td c plot cannot be associated with a point on the theoretical 
curve, because the observed data point cannot even be placed in this plot without knowing 
the 13 CO J = 1 — > line strength. Consequently, the correction factor in equation (40) is 
assumed to be unity (i.e., N d (H) = N dl (H) is assumed) for those positions outside the area 
mapped in 13 CO J = 1 — > 0. The latter problem is that the 13 CO observations alone do not 
yield accurate total masses for these clouds; the mass estimates from 13 CO J = 1 — > must 
be corrected for the missing coverage towards the edges of the clouds. These corrections 
were estimated by determining the total of the gas masses of the three big fields from the 
continuum observations for only those positions inside the I( 13 CO) = contour for each of 
the four cases (where Cases 1 and 2 are identical for the continuum-derived masses). These 
were then compared with the continuum- derived masses for all the gas in those fields (i.e., 
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the numbers listed as "Total" in Table 5). 

Before discussing this comparison, assessing the accuracy of this method is possible by 
comparing the continuum- derived masses inside the I( 13 CO) = contour with the 13 CO- 
derived masses in Table 4. These latter masses represent, in practice, the molecular gas 
mass inside the I( 13 CO) = contour plus the HI mass throughout the entire Orion fields. 
That is, the masses listed in the "Case 1", "Case 2", "Case 3", and "Case 4" columns of 
Table 4 are the 13 CO-derived molecular gas masses within the I( 13 CO) = contour plus 
the H I mass both inside and outside this contour. Therefore, to assess the accuracy of this 
correction for the lack of 13 CO coverage, the continuum-derived masses inside this contour 
must be corrected upwards for the H I mass that lies outside this contour. This H I mass is 
34000 M . After adding this mass to the continuum-derived masses inside the I( 13 CO) = 
contour, the continuum- derived masses agree with the 13 CO-derived masses of Table 4 to 
within 5 to 10%, depending on the case. 

Now to correct the total masses in Table 4 for the gas mass that is missing outside 
the I( 13 CO) = contour, we compare the continuum- derived masses within this contour 
with the masses in Table 5. This gives mass differences of 88000 M Q in the one-component 
cases and 87000 M in the two-component cases. Allowing again for the 34000 M of H I 
outside the I( 13 CO) = contour that is already implicit in Table 4 requires the masses in 
Table 4 to be corrected upwards by 54000 M in the one-component cases and by 53000 M 
in the two-component cases. The resultant best estimates for the Orion gas masses in the 
one-component case (for Case 2) and in the two-component case (for Case 4) are in Table 6, 
where the continuum-derived and 13 CO-derived masses agree to within about 5%. 

Table 7 shows that roughly 40 to 50% of the two-component gas mass is in the form 
of cold (T K < 10 K) gas (compare with Table 5). Consequently, the mass of gas or dust 
being detected at wavelengths much more sensitive to warm, rather than cold, gas or dust 
(e.g., in the Wien limit) must be corrected upward by a factor of at least 1.6. Figure 28 
shows that the cold component is on one edge of each of the small 2° HH-region fields (i.e. 
Orion Nebula and NGC 2024 Fields) and is spread out along the northern edge of the Orion B 
Field. There is also some cold dust towards the north-eastern edge of the A Orionis Field. 
A natural question at this point is about a realistic lower limit to the temperature: are 
gas and dust temperatures as low as 3 K really necessary for explaining the data? To test 
this, two-component models were applied, where T K1 was restricted to be > 10 K. This gave 
AT = — 7K and xl — 10.5. In other words, the minimum T<ji was still found to be 3K. 
Also, the x 2 is a factor of 1.9 higher than the best fit models that had no such restriction on 
T K1 . This is a worse fit at the 99% confidence level. Restricting both T dl and T K1 would have 
given an even worse fit. A less extreme restriction on T K1 of T K1 > 5K gives AT = — 2 K 
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and xl — 7-7. This fit is worse at almost the 90% confidence level and again the minimum 
T dl is 3 K. In short, avoiding having at least some dust or gas with temperatures near 3K 
is very difficult, provided that the basic assumption is valid. Further discussion of this is 
deferred to Paper III (Wall 2006a). 

Since the failings of the basic assumption could be relevant here, longer-wavelength 
observations confirming the existence of such cold gas or dust would be helpful. The longer 
wavelength FIRAS data were used to test this, but the results were inconclusive. The FIRAS 
data extend to a wavelength of about 4 mm, making them much less sensitive to the warm 
dust and, as such, much less insensitive to the cold dust. However, the FIRAS data also 
have a spatial resolution of 7°; the FIRAS beam solid angle is a factor of 100 larger than 
that of the DIRBE beam. Consequently, a direct comparison with the DIRBE data and the 
Orion fields as defined here is difficult. Nevertheless, such a comparison was attempted as 
follows: 

1. The spectra for 183 positions towards and around the Orion fields were obtained. Each 
position had two spectra: a low-frequency spectrum and a high-frequency spectrum. 
These were joined to form a single spectrum. Data in the channels that overlapped 
between the two spectra were averaged with l/cr 2 weights. This yields spectra that 
cover frequencies of from 68 to 2911 GHz (i.e. wavelengths from 4410 to 103 fim) with 
channels 13.6 GHz wide. 

2. The 60 highest frequency channels were then binned into 5 channels to improve the 
signal-to-noise ratio. The 12 channels that went into each of the 5 new channels were 
averaged with 1/a 2 weights. 

3. A background with the form Ibg = do + a\ csc(6 e //) was fitted to each channel from 
those positions outside the Orion fields and farther than half a FIRAS beam from 
the Galactic plane. The Oq represents an isotropic background component and is 
the cosmic background radiation. The a\ represents a scaling factor for the Galactic 
plane emission. The b e ff is the effective Galactic latitude for emission that follows a 
cosecant (b) law averaged within the FIRAS beam. Specifically, b e ff = arccsc((csc(b))) 
where (csc(b)) is the average of csc(b) within the 7° tophat beam of FIRAS. A good 
empirical approximation to b e ff is b e ff = \b\ — 5|6|~ L15 for |6| in the range of 4° to 
27°. After determining the best fit and a\ values for all the channels, the resultant 
background spectrum was subtracted from the FIRAS spectra centered on the Orion 
field centers. 

4. Two-component models of modified blackbodies were fitted to the 
background-subtracted spectra. These were of the form I v = r a (u / V ) /5 ] B u (T da ) + 
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Tbiy /vq) 13 B v (Tdb), where T a and T5 are the optical depths of components a and b, re- 
spectively, at reference frequency u Q . 

The best fitting curves had a warm component with temperatures between 17 and 20 K. 
The cold component was very hard to define. Assuming that this component even exists, 
the temperature is around 5 K, but could not be pinned down. Because of this temperature 
uncertainty, placing limits on the cold component mass was difficult. The cold-component 
mass was anywhere between and 10 times the warm-component mass. In any event, 
confirming the existence of this cold dust must await long-wavelength spectra with greater 
spatial resolution and sufficient spatial and wavelength coverage. 

4. Summary and Discussion 

We used the 240 /im continuum to the 13 CO J = 1 — > intensity ratio plotted against 
the 140 /im/240/im dust color as a diagnostic of the dust and gas physical conditions in 
the Orion A and B molecular clouds. Model curves were fitted to the data with orthogonal 
regression. Two sets of one-component models were applied to the data: Local Thermody- 
namic Equilibrium (LTE) and Large Velcocity Gradient (LVG). Two sets of two-component 
models were also applied to the data: two-component, LVG models and two-component, 
two-subsample, LVG models. In each case, the model fit gave the conversion of the observed 
intensities into gas/dust column densities and beam-area filling factors. The results of the 
fitting are summarized below: 

• One- Component, LTE Models: These models fit the data only poorly with a reduced 
chi-square of xl — 16.5. When considering systematic uncertainties, the only parame- 
ter to be fitted, AT, was very poorly determined: AT lies anywhere between —24 K and 
+6 K. Comparison of the continuum-derived gas column densities with those derived 
from the 13 CO J = 1 — > spectral line showed considerable scatter with xl = 16.6. 
The filling factors exhibit an understandable increase with increasing column density 
and are all less than unity. At best, the one-component, LTE models only crudely 
reproduce the overall trend in the data. 

• One-Component, LVG Models: These models also fit the data only poorly with a re- 
duced chi-square of xl — 16.9. When considering systematic uncertainties, the AT was 
again very poorly determined: AT lies anywhere between — 14 K and —1 K. Compari- 
son of the continuum-derived versus spectral line-derived gas column densities showed 
comparable scatter to the LTE model case, with xl = 15.5. The filling factors behaved 
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similarl to those in the case of the LTE modeling. The LVG model results overall are, 
at best, only marginally better than the single-component, LTE models. 

• Two- Component, LVG Models: In contrast to the one-component models, these models 
explain the overall trend in the data with xl — 5.7; this is better than the one- 
component models at a confidence level of better than 99.9% according to the F- 
test. Also in contrast to the one-component model results, AT is known to be within 
1 K of K even when the systematic uncertainties are considered. The component-0 
dust temperature seems to be accurately known and is consistent with the large-scale 
interstellar radiation field, i.e. 18 K. Other parameters, such as gas densities and 
column densities per velocity interval, are known from factors of a few to more than 
an order of magnitude about some geometric mean. The continuum-derived versus 
spectral line-derived gas column densities showed much better agreement for these 
models than for the one-component component models, the former having xl — 7.4. 
Despite their advantages, these models suffer from two problems. The first is that 
many of the filling factors have unphysical values — i.e., greater than unity. The second 
problem is that the model curves do not fit well to the data with 140 /im/240 /im color 
temperatures greater than 20 K. 

• Two- Component, Two- Sub sample, LVG Models: Here the sample of data points was 
split into two subsamples: one with dust color temperatures < 20 K and the other with 
such temperatures > 20 K. These models have all the advantages of the simpler two- 
component models described previously, and also fix the problems of the unphysical 
filling factors and the poor fit to the data points with 140//m/240/xm color temper- 
atures > 20 K. The only disadvantage is that this approach is ad hoc. Nevertheless, 
the derived parameter values agree with those of the simpler two-component models 
to within the systematic uncertainties. 

• Cloud Masses: After correcting for systematic effects, such as the lack of coverage of 
the 13 CO J = 1 — > map, the total hydrogen gas mass derived in the Orion fields is 
4.2 x 1O 5 M , as found in the more realistic two-component case (and is 2.5 x 1O 5 M 
for the less realistic one-component case). Due to the systematics, the uncertainty in 
this mass is ±40%. The extra mass in the two-component case is due to the apparent 
presence of cold dust/gas with temperatures of 3-10 K. Attempting to confirm the 
existence of this cold dust/gas with the long- wavelength FIRAS data was inconclusive 
due to the insufficient spatial resolution. 

The above results have a number of astronomical implications, which are mentioned 
here only briefly: 
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1. Dust Temperatures Equal to Gas Temperatures: If this surprising result applies to the 
gas and dust on large scales in our Galaxy, then models of gas-to-dust thermal coupling 
must be revised. Typically, these models assume spherical grains with a single large 
grain size and a smooth grain surface (e.g., see Goldsmith 2001; Burke & Hollenbach 
1983). Clearly, relaxing these simple assumptions will increase the surface area-to- 
volume ratio of the grains and improve the dust-gas thermal coupling, reducing AT 
(see Papers III and Ilia for more details). Having dust and gas temperatures equal on 
galactic scales implies that large-scale molecular gas temperatures are nearly double 
the temperatures previously believed (i.e. closer to 20 K than to 10 K). Also, one dust 
grain alignment mechanism would be ruled out. 

2. At least two temperatures are needed along each sightline to give reasonable column 
densities: The two-component models show good agreement between spectral line- 
derived and continuum-derived column densities. The main difference between the 
two components is that one component has unvarying temperature from sightline to 
sightline and the other component has varying temperature. Therefore, except for those 
few sightlines where the two components have the same temperature, each sightline 
has two temperatures (i.e. one temperature — gas and dust — for each component). 
Requiring more than one temperature along each sightline is supported by the work of 
Schnee et al. (2006), who find that assuming a single temperature along each sightline 
does not give reliable column densities in the Perseus and Ophiuchus molecular clouds. 

3. Cold Gas/Dust: If similar two-component models are correct and are applicable to 
other large clouds in the Galaxy, then the total mass of molecular gas in the Galaxy 
must be revised upwards by about 60%. 

4. Millimeter Continuum to 13 CO J = 1 — > Line Ratio as a Temperature Diagnostic: 
For the millimeter continuum to 13 CO J = 1 — > intensity ratio above some threshold, 
this ratio can serve as a diagnostic of the gas/dust temperature, even in the high- 
temperature limit. 

5. X-Factor: An improved explanation for the X-factor is applicable. This will be dis- 
cussed in detail in Wall (2006b). 

Discussing the above results in depth would be premature at this point. We must first 
examine the modeling and various systematic effects more carefully. Even though we have 
tested the effects of some systematics on the model curves (i.e., by applying scale factors 
or changing the starting grid when fitting the models), we need to find and understand any 
inherent biases or flaws that the model fitting may have. For example, the two-component 
models consistently find a AT equal to or near K, in spite of the scale factor applied to 
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the data. Is this because the dust-gas temperature difference really is zero? Or is it simply 
because these models are biased towards AT = K despite its true value? And what of 
the reliability of the other parameter values? The component-1 density, for example, is 
supposed to be known to within a factor of a few of 600 i/ 2 cm~ 3 , according to the simple 
two-component models (see Figure 20). And yet, the two-subsample, two-component models 
suggest densities that are one to three orders of magnitude higher than this. Clearly, the 
methods used thus far to examine the systematics are incomplete. Consequently, Paper II 
examines the results of running the models on simulated data and will answer the concerns 
raised here. 

In addition to those concerns, there are other questions not yet addressed by the methods 
used here. Specifically, what is the effect of the background subtractions used? How will 
dust associated with HI affect the results? Does changing the spectral emissivity index f3 
appreciably affect the results? Are there alternative kinds of models that would also explain 
the data? The cloud positions modeled here only represent 26% of the area of the Orion 
clouds; how representative are the results of the clouds as a whole? Paper III examines these 
questions. 

All these questions and concerns are addressed in Papers II and III. Paper III also 
discusses the scientific implications of the results in more detail. 
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A. Column Densities from Optically Thin Rotational Lines of CO 

In this appendix, the relationship between the molecular gas column density and inte- 
grated line strength of an optically thin rotational line of CO, or one of its isotopologues, is 
derived. The derivation detailed below is more appropriately applied to the isotopologues 
of CO, such as 13 CO and C 18 0, than to CO itself (i.e. 12 C 16 0), because the low-J lines 
of CO (i.e. J = 1 — > 0, 2— >1, 3^2) are usually optically thick, whereas the correspond- 
ing lines of its isotopologues are usually optically thin for molecular gas on the scales of 
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many parsecs. Such a derivation is helpful for a number of reasons. One reason is that this 
relationship is central to this paper and knowing which approximations were not used, as 
well as which ones were used, is important. Another reason is that the full derivation is 
difficult to find in the literature. Furthermore, aside from the usual assumptions made in 
determining column densities (i.e., optically thin line, non- varying gas physical conditions 
along the line of sight, and sometimes LTE), the equations in the literature are sometimes 
based on additional invalid or unnecessary assumptions. For example, the Rayleigh- Jeans 
approximation is sometimes used and is both invalid and unnecessary. (Ironically enough, 
using the Rayleigh- Jeans approximation actually results in a more complicated expression 
than not using this approximation.) 

The optical depth of the J — > J — 1 rotational transition, r^(J — > J — 1), is related to 
the volume absorption coefficient, «„(J — > J — 1), by 

r^J — J-l) = Jk„(J^ J-l) ds (Al) 

or 

r„(J - J - 1) = ^- J (B a _ M n^ - Bj t n^hv^MJ - J - 1) ds . (A2) 

The integral is over the line of sight through the cloud. The n J _ 1 and n, are the volume 
densities of those CO molecules in the J— 1 and J rotational levels, respectively. The B J U 
and B a J _ 1 are the Einstein B-coefficients defined in terms of the mean intensity of the 
radiation field, rather than in terms of the energy density. The Bj_i j is the absorption 
rate coefficient for absorptions from level J— 1 to level J. The Bj 5 j_i is the rate coefficient for 
stimulated emission from level J to level J— 1 and is the stimulated emission correction to the 
optical depth in equation (A2). The hi/ J31 is the energy difference between levels J— 1 and 
J, where h is Planck's constant and v } T _ 1 is the frequency of the J — ^ J — 1 transition. The 
quantity </>„(J — > J — 1) is the line profile function in terms of frequency, v. It is normalized 
such that 

J <f> v (J -> J - 1) du = 1 . (A3) 

spectral line 

The integration in (A2) converts the volume densities, n 7 _ 1 and n a , into the corresponding 
column densities, N J _ 1 and N j; yielding 

r„(J - J - 1) = ifB^N,., - B^N^h^.^J -J-l) . (A4) 

Note that the assumption of LTE has not been used. Similarly, assuming constant gas 
physical properties along the line of sight has not been applied. (Strictly speaking, the line 
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profile function can change along the line of sight. Nevertheless, the 4> U (J — > J — 1) that 
appears in A4 is the appropriately weighted average along the line of sight of the 0„(J — > J— 1) 
that appears in A2.) The following relations eliminate the Einstein B-coefficients in favor of 
the spontaneous emission rate coefficient, A 7 7 _ 1 : 

= — B,,,., (A5a) 



and 



c 2 



= A w-i > ( A5b ) 



2hz/ 3 

J.J-l 



where c is the speed of light and g T and g J _ 1 are the degeneracies of levels J and J— 1, 
respectively. Using (A5a) and (A5b), equation (A4) becomes, 



r„(J - J - 1) = — — gA,-! ( — " — J ^ J " 1) • (A6) 



The frequency profile function, <f) u (J — > J — 1), is replaced with the velocity profile function, 
V (J — > J — 1). Since (f) u (J — > J — 1) df = V (J — > J — 1) dv, 0^(J — > J — 1) = A, J _ 1 V (J — ► 
J — 1), where A T J _ 1 is the wavelength of the J — > J — 1 transition. Substituting this into (A6) 
yields 

r v (J -> J - 1) = " |A " ;A:i ^ fl - i^g£^ v (J _ J _ 1) . (A7) 
The excitation temperature between levels J and J— 1, T X (J — > J — 1), is defined by 

= exp -— — — , (A8) 



kT x (J -> J - 1) 

where k is Boltzmann's constant. Up to now, the assumption of constant gas physical 
properties along the line of sight and the assumption of LTE have not been used. 

The term "excitation temperature" is often used incorrectly in the literature. Under 
the assumptions of some optical depth limit and LTE, a ratio of two molecular lines gives 
a temperature. In the LTE limit, this line-ratio temperature is the kinetic temperature 
of the gas. However, it is often called the "excitation temperature of the g though 
excitation temperature is a physical property of the gas itself, like kinetic temperature, 
rather than a parameter determined by the physical properties. The excitation temperature 
properly represents the population ratio of two particular levels of a particular molecule 
and not a physical property that characterizes all the levels of all the molecules in a region 
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of a molecular cloud. A line-ratio temperature represents three or four levels and is not 
an excitation temperature, but depends on the two excitation temperatures of the lines 
forming the ratio (as well as their optical depths). In general, temperatures from line ratios 
and the excitation temperatures of the various transitions are numerically distinct from 
one another, except in the limits of either low or high density. In the former case, all 
excitation temperatures equalize to the temperature of the ambient radiation field (i.e., 
the cosmic background temperature in the absence of other sources) and the temperatures 
derived from line ratios are equal to this temperature. In the latter case, all the excitation 
and line ratio temperatures are the same as the kinetic temperature of the gas. As such, 
many articles actually mean "line-ratio temperature" or "kinetic temperature" when they say 
"excitation temperature" . In particular, the term "kinetic temperature" is appropriate when 
such articles compute column densities from a single optically thin rotational line of CO. 
Unless a non-LTE radiative transfer code (e.g., the Large Velocity Gradient or LVG code) 
is used, along with some constraints on the gas physical conditions, to explicitly consider 
the excitation temperature variation with position in the rotational ladder, the excitation 
temperatures are implicitly assumed to be numerically constant throughout the rotational 
ladder in computing the partition function. Therefore, this adopted and unvarying (with J) 
excitation temperature is actually the adopted kinetic temperature (or the temperature of 
the ambient radiation field in the case of very low density). 

Now the gas physical conditions are assumed constant along the line of sight, so that 
N J /N J _ 1 = n J /n J _ 1 . Substituting (A8) into (A7), integrating over the line velocity profile, 
and rearranging gives 



N, 



87T gj 



A 3 



1 — exp 



hz/ 



j, j-i 



kT x (J -J-1) 



-i 



J r v (J - J - 1) dv , 



(A9) 



spectral line 



where we used the normalization condition (A3), but applied to V instead of <f> v . Following 
Kutner and Ulich (1981), the radiation temperature, T R , which is the brightness temperature 
defined in the Rayleigh- Jeans sense (as opposed to the Planck sense), is given by 



T R (J-J-1) = 



hi/, 



exp 



hz/ 7 



kT x (J -J-1) 



- 1 



exp 



A", 



kT 



BG 



-i -1 



1 — exp ( — r v ( J — J — 1) ) , (AlO) 



where T BG is the cosmic background temperature. Note that, even though this is a Rayleigh- 
Jean brightness temperature, the Rayleigh- Jeans approximation is not used. For optically 
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thin lines, (A10) can be rewritten as 



kT R (J-> J-l) 
hi/, , , 



cxp 



hz/ T 



kT x (j -> j - r 



- 1 



r v (J -> J - 1) C E 



(AH) 



where C BG is a correction for the cosmic background radiation and is given by 



exp 



(■ 



hi/. 



a 



i - 



\kT x (J^J-l) 



- 1 



cxp 



kT 



BG 



(A12) 



- 1 



Notice that this correction factor is zero when T X (J — > J — 1) = T BG , because a line with 
an excitation temperature equal to that of the cosmic background is not seen against this 
background. 

Many expressions in the literature for CO column density implicitly assume that C BG is 
unity. This is a reasonable assumption when T X (J — > J — 1) is appreciably greater than T BG . 
For example, for the 13 CO J = 1 -> line (u w = 110.201 GHz) the T x (l -> 0) is probably 
close to the kinetic temperature, T K , which is often no lower than 10 K in most cases. For 
T K = 10 K, C BG = 0.88, thereby affecting the column density derivation by about 12%. 
For most purposes, this accuracy is sufficient and is even better at higher temperatures: 



for T K = 15 K, C 



bg 



0.93 and for T K = 20 K, C 



BG 



0.95. The accuracy is also better 
at higher frequencies: for the 13 CO J = 2 — > 1 line (u w = 220.399 GHz), again adopting 
T x (2 -> 1) = T K for simplicity, the C BG = 0.96, 0.98, 0.99 for T K = 10, 15, and 20 K, 
respectively. Nevertheless, for the work described in this paper, the C BG correction was 
included in the computations. 

Rearranging (All) gives 



1 — exp 



hz/ T 



kT x (J -> J - 1) 



Substituting (A13) into (A9) results in 



N, 



! ' .t-l\j-l nZ/ .T,.T-l CbG 



exp 



-l 

r v (J -> J - 1) 
kT R (J^J-l) 



ni/ J,J-l ^BG 



hz/ T 



cxp 



kT x (J 



1) 



hu 



kT x (J -> J - 1) 



(A13) 



) / T R (J-J-1) dv 



(A14) 
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The excitation temperature, T X (J — > J — 1), is not known in general, but is now eliminated 
by using (A8) and the expression N J /N J _ 1 = njn } _ 1 : 



8tt k f 
Nj = /T R (J^J-l)dv . 

Vl A ?,J-l hZ/ J,J-l C BG J 



(A15) 



Even though T X (J — > J — 1) is still implicitly present in C BG , the dependence of C BG on 
T X (J -> J - 1) is very weak for T X (J -> J - 1) » T BG . 

Equation (A15) can be simplified further using 

64 * 4 "?..,-, I J A 2 



A »- = ^UnJ" • (A16) 

where \x is the permanent dipole moment of the molecule. This is 0.11 Debye for CO, where 
1 Debye = 10- 18 esu-cm in CGS units or 3.33564 x 10" 30 m-C in SI units (see CRC Handbook 
of Chemistry and Physics 1998-99). The energy above the ground state of level J, hz/ J0 , is 
well represented by 

hz/ J0 = hBJ(J + l) , (A17) 

where B is the rotational constant of the molecule (not to be confused with the Einstein 
B-coemcients). From (A17) it follows that 

z/. M _ 1= 2BJ . (A18) 

Equations (A17) and (A18) assume that the molecule is a rigid rotor and neglect the effects 
of centrifugal stretching. Nevertheless, (A18) gives the correct observed frequencies of the 
rotational lines of CO up to the J = 10 — > 9 line to within 0.06% or better. (Note that this 
is still very inaccurate for obtaining radial velocities during actual observations, but is more 
than sufficient accuracy for column density determinations.) Using (A18), the rotational 
constant can be easily determined from B = v 10 /2, where v w is the observed frequency of 
the J = 1 — > line. For CO and its isotopologues, these values are B(CO) = 57.64 GHz, 
B( 13 CO) = 55.10 GHz, and B(C 18 0) = 54.89 GHz. Equation(A16) is now written as 

512tt 4 B 3 f J 4 \ 2 



A = /i 2 . (A19) 

3hc 3 V2J + 1/ 



Substituting equations (A19) and (A18) into (A15) yields the useful expression 

3k(2J + 1) 

Nj = — Ij,i , (A20) 

J 16vr 3 BJ 2 /i 2 C BG 
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where I, J _ 1 = J T R (J — > J — 1) dv is often called the integrated intensity of the line (although 
it is really the velocity- integrated radiation temperature of the line). Normally, CGS units 
are used, but I T 7 _ 1 is in units of K • km • s _1 ; consequently (A20) becomes 

3k(2J + 1) 

Nj — — ; -e— ' • lO'-l;.; iK - km-s- 1 ) , (A20a) 



16tt 3 BJ V C 



BG 



where all quantities with units, except I,,^, are in CGS units. So far, the gas physical 
conditions have been assumed constant along the line of sight. LTE has not yet been assumed. 
Therefore, for general non-LTE conditions, the total column density N(CO) can be estimated 
by observing a number of optically thin rotational lines of CO, using (A20a) to determine 
the column density in the upper levels of the observed transitions, and sum the Nj values. 
This sum must then be corrected for the column densities in the upper levels of unobserved 
transitions and for the column density in the ground state, N . (Notice that A20a is not 
defined for N .) As such, the fraction of the total column density in the ground state and in 
the upper levels of unobserved lines must be estimated by modeling the physical conditions 
from the observed line ratios. This fraction then corrects the sum of the observed N a values 
to an estimate of the true N(CO). 

When only one or two optically thin lines of CO are observed, LTE is often assumed. 
Notice from (A19) that the spontaneous emission rates are lowest for the lowest-J transitions, 
implying that these lines are closest to LTE. So if the lines observed are low-J lines like 
J = 1 — > or J = 2 — > 1, then the LTE assumption might reasonably estimate the total 
CO column density from the integrated intensity of a single observed line. The total column 
density of CO is given by 

oo 



N(CO) = = N, ( -M £ / • ( A21 ) 

J'=0 ^ iN| .j/ J'-n iN| 



7^ N 



The J in this expression is the upper-J of the observed transition. Replacing all occurrences 
of J— 1 in (A8) by and, because LTE is now assumed, replacing T X (J — > J — 1) by the 
kinetic temperature, T K , for all J, converts (A21) to 

N(C0) = 2lTT Q(TK)exp (kf J ' (A22) 

(A22) uses g, = 2 J + 1 and g = 1. The Q(T K ) is the rotational partition function given by 

* , / hBJ'(J'+i)\ 
Q(T K ) = ^(2J' + l)exp A- L) , (A23) 

J'=0 ^ kl K / 



-52 - 



where (A17) was used. When the kinetic temperature is high enough, the summation vari- 
able, J', can be treated as a continuous variable (i.e. real instead of integral) and the 
summation in (A23) is treated as an integral over J'. This gives the simplified, but approxi- 
mate, 

kT 

Q(T K ) — ■ (A24) 

hB 

For CO, 13 CO, and C 18 0, (A24) is T K /2.766, T K /2.644, and T K /2.634, respectively. 

As is the case for the cosmic background correction factor, C BG , the approximate form 
of the partition function, (A24), is often used in the literature. Again, this approximate 
form is appropriate most of the time, because usually T K >> hB/k. For T K = 10, 15, and 
20 K, the ratios of the approximate Q(T K ), given by (A24), to the exact Q(T K ), given by 
(A23), are 0.91, 0.94, and 0.96, respectively. Again, for the work described in this paper, the 
approximate form was not used. It is easy to see why: the effects of neglecting the cosmic 
background correction and using the approximate rotational partition function reinforce each 
other and result in a theoretical I( 13 CO) (for the J = 1 — > line) that is too high by 14% at 
T K = 15 K. This is non- negligible for the work here. 

Substituting (A17) and (A20) into (A22) results in 

3k /hBJ(J + l)\ 

N(CO) ^6,3B JV C Ba Q(T - )eXP (^H) ■ (A25) 

Using the same units as in (A20a) gives 

3 x 10 5 k /hBJ(J + l)\ 

n(co) = i WCm q(T « )exp (^H K ■ km ■ s ) • (A25a) 

where N(CO) is in units of CO molecules per cm 2 . Note that the C BG that appears in 
(A25) and (A25a) is now a function of T K instead of T X (J — > J — 1). Substituting in the 
fundamental physical constants gives 



Q(T K ) / BJ(J + 1)\ 
exp 4.80 x 10 I,, , K' km-s 



N(CO) = 6.90 x 10 24 exp ( 4.80 x lO' 11 — ' - ) I T t-1 ( K ■ km ■ s' 1 ) . (A26) 

BJ 2 C BG V T K / 

Choosing a particular isotopologue of CO, 13 CO, then applying the value for the rotational 
constant for 13 CO results in 



N( 13 CO) = 1.25 x 10 14 ^ exp ( 2.64 ^±H) I^^CO) , 

J ^BG V *-K / 



(A27) 
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where Ij J _ 1 ( 13 C0) is still in units of K- km-s -1 . Equation (A27) is used in the current paper 
for determining the molecular gas column densities in the LTE case. When T K >> 2.644 K, 
(A27) can be simplified using the asymptotic form for Q(T K ), (A24), and C BG ~ 1 so that 

N( 13 CO) = 4.73 x 10 13 ^ exp^2.64 J ^ T +1 ^ Ij^^CO) . (A28) 

Note that (A28) was not used in the current paper, but is still a useful approximation: the 
upward corrections for this equation are 14%, 6.5%, and 1.8% for T K = 15, 30, and 100 K, 
respectively. 

Using the appropriate abundances with (A26) and (A27) allows determination of N(H 2 ). 
Reasonable values for these abundances are X( 12 CO) = 8 x 10~ 5 , X( 12 CO)/X( 13 CO) = 25 to 
65, and X( 13 CO)/X(C 18 0) = 5 to 6 (e.g., Encrenaz et al. 1975; Dickman 1975; Knapp & Jura 
1976; Frerking et al. 1982; Langer k Penzias 1990; Taylor & Dickman 1989). Note that the 
abundances of the isotopologues near cloud surfaces can be different from the numbers given 
here, which are the abundances for cloud interiors (see, e.g., Taylor & Dickman 1989; Frerking 
et al. 1982). Simply dividing the CO column density by the abundance can give a reasonable 
estimate of N(H 2 ), but Frerking et al. (1982) find that a zero-point N(H 2 ) must also be added 
to explain their observations of the Taurus and p Ophiuci molecular clouds. Such a zero-point 
column density was not used in the current work because it varies by a factor of a few from 
cloud to cloud (see Frerking et al. 1982) and because it would make little difference to the 
results and conclusions. For the 13 CO abundance, we adopt X( 12 CO)/X( 13 CO) = 63. This 
is an average of the value observed at the position (0.5, 2.5) observed by Langer & Penzias 
(1990) and the interpolated value expected for the galactocentric radius of the Orion clouds 
(see their Figure 6). This average is probably more representative of the Orion clouds as a 
whole and only differs by 6% from the observed value at position (0.5, 2.5). Note that the 
observed value at the (0,0) is probably peculiar because it lies within a region with strong 
UV fields (see Langer & Penzias 1990). 

Briefly discussing other approaches for determining column densities with optically thin 
CO lines is worthwhile. For example, sometimes the Rayleigh- Jeans approximation is used. 
In that case, the excitation temperature of the transition is assumed to be sufficiently large 
that the approximation hv } J1 /kT x (J — > J — 1) << 1 is valid. Equation (All) would have 
the much simpler form, 

T R (J -> J - 1) ~ T X (J -> J - 1)t„(J -> J - 1) (A29) 

However, adopting this simpler expression actually complicates the final formula (i.e., A25) 
for the column density because the [1 — exp(— \w } J :1 /kT x (J — > J — 1))] _1 factor is not 
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replaced with a simpler form. Consequently, the LTE expression (A25) has an additional 
factor, (hz/j j _ 1 /kT K ) [exp(hz/j J1 /kT K ) — l] -1 . For some reason, this additional factor is 
retained, even though the Rayleigh- Jeans approximation used to validate the use of (A29) 
requires its replacement with unity. Because the high-T K limit applies, Cbg is set to 1 and 
the approximation for the partition function, (A24), is used. The latter allows cancellation 
of the extra T R : . The resultant expression reasonably estimates the column density only in 
the high-T K , low- J limit. The ratios of the column density derived from the expression using 
the Rayleigh- Jeans approximation to that of (A25) for the 13 CO J = 1 — > line are 0.61, 
0.79, 0.86, and 0.91 for T K = 10, 20, 30, and 50 K, respectively. For the 13 CO J = 2 -> 1 line 
the corresponding ratios are 0.49, 0.72, 0.80, and 0.88. For the 13 CO J = 3 -> 2 line these 
ratios are 0.37, 0.63, 0.74, and 0.83. If column density estimates within a factor of 2 of those 
from (A25) are considered acceptable, then, even for the rather common kinetic temperature 
of 10 K, the estimates based on the Rayleigh- Jeans approximation are unacceptable for all 
transitions higher than J = 1 — > 0. Even for J = 1 — > 0, the column density estimate 
must be corrected upwards by 64%. A similarly large upward correction, 59%, is necessary 
even for a kinetic temperature of 20 K for the J = 3 — > 2 line. Avoiding the Rayleigh- Jeans 
approximation keeps the final expression more accurate and less complicated. (Of course, 
if the Rayleigh- Jeans approximation is applied correctly, then the additional factor is unity 
anyway. The resultant expression would be correct in the high-T K limit, where C BG is set to 
unity, which is appropriate for this limit.) 

Another approach is to use the corresponding line of 12 CO, which is optically thick, 
assume LTE and that the gas fills the beam in the velocity interval about the 12 CO line 
peak, and solve for the optical depth r( 13 CO) or r(C 18 0) from the 13 CO/ 12 CO or C 18 0/ 12 CO 
intensity ratios. In addition, the 12 CO line peak intensity estimates the gas kinetic temper- 
ature. ( Note that this is often inappropriately referred to as the "excitation temperature" . 
It is inappropriate because LTE was assumed in computing r and in the expression for the 
column density.) The optical depth, kinetic temperature, and line velocity-width are then 
used in an LTE expression for the column density that assumes the high-T K limit for the 
rotational partition function. This method has the advantages that assuming the optically 
thin limit is not necessary and that the additional information provided by the 12 CO line is 
used. However, it also has the disadvantages that the high-T K limit for the partition func- 
tion is used and that the gas fills beam in the interval about the 12 CO line peak is assumed. 
The former assumption does not produce very large errors: even down to T K = 2.8 K, the 
upward correction to the column density found from using the full partition function is only 
40%. The latter assumption is likely wrong for external galaxy observations, because the 
molecular gas is unlikely to fill the beam in any single velocity interval. (As discussed in 
this paper, it is also likely wrong for observations of Galactic molecular clouds, but may 
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still be roughly correct to within a factor of about 2.) In fact, this assumption can lead, in 
certain cases, to column density estimates that are wrong by factors of a few or more. If we 
consider a few cases with I(C 18 J = 1 — > 0) a factor of 20 down from I( 12 CO J = 1 — > 0) 
and if T R ( 12 CO J = 1 — > 0) is weak, then the column density estimates are too high by 
factors of 2 to 5. For example, if the filling factor in each velocity interval is 0.1 and if 
T R ( 12 CO J = 1 — > 0) = 0.1 K, the column density estimate is too high by a factor of 5. 
Even in the less extreme case of T R ( 12 CO J = 1 — > 0) = 0.5 K and a filling factor of 0.2, the 
estimate is too high by a factor of 2. Even if T R ( 12 CO J = 1 — > 0) is much stronger, 10 K, 
the estimate is still a factor of 2 too high for a filling factor of 0.1. Despite this shortcoming, 
the correction factor to the column density estimate is often much less extreme than the 
filling factor itself. In any event, this method is easily improved by appropriate application 
of the filling factor to the equations used to derive the optical depth, kinetic temperature, 
and the column density itself. 

One final point is that elaborate non-LTE radiative transfer models suggest that LTE- 
derived CO column densities can be underestimates of the true column densities by factors 
of 1.3 to 7 (Padoan et al. 2000). As shown in the current paper, the assumption of LTE 
can indeed underestimate the true column density; similar factors are derived using much 
simpler two-component models that employ the LVG code. 



B. Ratio of the Submillimeter or Millimeter Continuum to an Optically Thin 

13 CO Rotational Line in the LTE Case 

The optically thin continuum emission of the dust has a specific intensity, L,, given by 

I„ = r„B„(T d ) . (Bl) 

The dust optical depth r v is given by 

t v = /c l/ x d // H m H N(H) , (B2) 

where k v is the dust mass absorption coefficient at frequency u, ji H is the correction for 
the presence of helium, x d is the dust-to-gas mass ratio (where "gas" means hydrogen and 
helium), m H is the hydrogen atom mass, and N(H) is the column density of hydrogen nuclei. 
Specifically, 

N(H) =N(HI) + 2N(H 2 ) , (B3) 

where N(HI) is the column density of atomic hydrogen and N(H 2 ) is the column density 
of molecular hydrogen. Substituting (B2) into (Bl), using the explicit form for the Planck 
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function, and converting from cgs units to MJy • sr 1 yields 

hz/ 3 /t. 



I„( MJy ■ sr' 1 ) = 2 x 10 17 ^! - N(H) . (B4) 



c 2 [exp(^) - 1] 



The mass absorption coefficient at frequency v expressed in terms of that coefficient at a 
reference frequency u is 

/ v \ P 

Kv= [ — ) Kvo , (B5) 

where (3 is the spectral emissivity index of the dust. Substituting (B5) into (B4) gives 

I,(MJysr- 1 ) = 2xl0 17 ^ H m H 

^[exp(^)-l] 

Now that the continuum intensity is expressed in terms of temperature and column 
density, the same must be done for the 13 CO rotational line. In the LTE, optically thin case, 
such an equation is simply the rearrangement of (A25a). At this point the assumption N(HI) 
= is adopted because the molecular gas column densities dominate for the majority of the 
positions studied here. (Nonetheless, there are positions for which the atomic gas column 
density is appreciable, introducing possible complications in the interpretion of the results. 
These complications are examined in detail in Paper III.) Applying the abundance of 13 CO 
to change N( 13 CO) to N(H 2 ), using (B3) with N(HI) = to change N(H 2 ) to N(H), and 
rearranging (A25a) results in 

n3 _ 87r 3 Bj 2 /i 2 c BG x( 13 co) ( hBj(j + r 



W»CO)= 3X — ° Q(TJ ' e^-^^j N(H) , (B7) 
where Ij 5 j_i( 13 CO) is in units of K • km • s _1 . Dividing (B6) by (B7) gives 

/ hBJ(J+l)\ 

I,(MJy ■ sr' 1 ) hkz/ 3 +V x d/ i H m H Q(T K ) ex V{^r) 

Vi Ij,J-i( 13 CO) • vr 3 BJV^C BG X(i 3 CO)c% xp ^)_ 1 ' 

where is in units of MJy • sr -1 • (K • km • s -1 ) -1 . Substituting in the fundamental physical 
constants and for the rotational constant of 13 CO, 

„ I 2.64J(J+1) 

= 1-67 x 10 4 S^Z^^ 7,, t~ \ '- ■ (B9) 



A^ (;Um )j 2 x(i 3 CO)C BG f i^ ).! 



where \ (/j,m) and X(/im) are the wavelengths in [im corresponding to frequencies z/ and u, 



respectively. Note that r A . is still in units of MJy-sr 1 -(K-km-s 1 ) 1 . Although not explicitly 
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stated in every case in W96, that paper adopted the following values: k u q = 40cm 2 /g for v 
corresponding to A = 100 ^m, x d = 0.0075, fi H = 1.4, and (5 = 2. X( 12 CO)/X( 13 CO) = 63 
was adopted for the 13 CO abundance (see Appendix A), so that equation (B9) becomes 

/ 2.64J(J+1)\ 

OCT ) ex P — t — I 
r Aj = 5.51 x W™ —^±d- / \ J . (B10) 

A5 W J2C Boexp(^)-l 
If A is set to 240 /im and J to 1, then 

OCT ) ex p(|^ 
r 240 = 69.2^ V \ K ^ ^ (Bll) 

C BG exp ^j_! 

where the second subscript on the r (i.e. the J) was dropped for clarity. Equation (Bll) is 
used in the current paper for the LTE modeling of the continuum to line ratio as a function 
of T d . 

It is interesting to examine the behavior of r 240 in the low- and high-temperature limits. 
For simplicity we assume that T d and T K are the same and equal to some T. We also assume 
that T >> 2.644 in both the low- and high-T limits, so that the simpler form of the partition 
function, (A24), can be used and so that C BG ~ 1. The low-T limit means that T << 59.9 K 
and so 



r 



/ 54 6\ 

26.2Texp(-^j . (B12) 

The upward correction to this formula is 16% and 23% at T=15K and 30 K, respectively. 
The high-T limit is where T >> 59.9 K and we have 

r 240 = 0.437 T 2 . (B13) 

The downward correction to this formula is 20% at T=100K. In short, the continuum to 
line ratio has a strong Wien-like rise at low temperatures and a relatively weaker T 2 rise 
at high temperatures. This high-temperature rise is still very strong compared to that of 
ratios of CO rotational lines from the same isotopologue and to that of ratios of continuum 
intensities at two different wavelengths. 
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C. Uncertainty Analysis of the LTE Derivation of the 13 CO Column Density 

To derive the uncertainty in the 13 CO column density in the LTE case, we start with 
equation (A27) with J=l: 

N( 13 CO) = N Q(T K ) C^(T K ) exp (J^j I( 13 CO) , (CI) 

where N = 1.25 x 10 14 and T 10 = hz/ 10 /k = 5.28 K. The uncertainty in the derived column 
density, 0"(N), is given by 



Substituting equation (CI) into (C2) and then dividing by (C2) yields 

<x 2 (N) /T K Q'(T K ) T K C' BG (T K ) T 1() \ 2 a\T K ) , a\\) 



N 2 V Q(T K ) C BG T K ; Tl + P • (C3) 

Differentiating (A12) gives C^ G (T K ): 

^ exp (^) 

C; G (T K )= ~ \«' . (C4) 
exp fe) - 1 

Similarly, differentiating (A23) gives Q'(T K ): 

q^ t k) i)(2j + i} 6x p(- BtJ t + 1} ) ' (cs) 

K J=l V / 

where B T = hB/k = 2.64 K and is the molecular rotational constant in units of Kelvins. To 
simplify the above expression, only the first two terms in the summation are necessary at 
low temperatures, i.e. T K < 2B T , to achieve 5% accuracy or better. At high temperatures, 
i.e. T K > 2B T , the approximate expression (A24) can be used. This results in the following 
approximations for Q'(T K ): 

B A ( B T J(J + 1)\ 
Q( T k) = ^^J(J + l)(2J + l)exp [-^- -), T K <2B T (C6a) 

K J = l V J-K / 

Q'(T K ) = J- , T K >2B T . (C6b) 
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The error in this approximation is no worse than 5% for T K just below 2B T or 5.28 K and 
better than this at higher and lower temperatures. The uncertainty in T K itself, cr(T K ), is 

a 2 (T K )=<7 2 (T d ) +( 7 2 (AT) . (C7) 

Therefore, the uncertainty in the derived 13 CO column density is given by equation (C3), 
with equations (C4), (C5), (C6b), (C6b), and (C7) substituted in. 



D. Uncertainty Analysis of the Column Density Derived from the Dust 

Continuum Emission 

The Wien limit is used for estimating the uncertainties in the dust-derived column 
densities. Combining equations (Bl) and (B2) and using the Wien approximation gives 

N = I„ D exp , (Dl) 



T. 



where D = [(2hz/ 3 /c 2 )/t iy/ u H m H x d ] 1 and F = hu/k. The uncertainty in the dust-derived 
column density, cr(N), is 



^ = {wj ^(jm) ff2(Td) (D2) 



Evaluating (D2) using (Dl) and then dividing by (Dl) results in 

^ 2 (N) _ <x 2 (U a 2 (T d ) F 2 
N 2 11 + T 2 T 2 • 1 ' 

For v corresponding to A = 240 /xm, F = 60.0 K. Equation (5) gives a(T d ). 



E. Uncertainty Analysis of the LVG Derivation of the 13 CO Column Density 

As stated in Section 3.2, the column density of molecular gas derived from the 13 CO J = 
1 — > line, Ni 3 (H 2 ), using LVG model results can be written as 



y = R(ai, a 2 , a 3 ; x 1 ,x 2 , x 3 ) x 3 



(El) 
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where 

y = N( 13 CO) (E2) 

and 

R(a a a-x x x ) ~ ^^/^od (m 

R(a u a 2 , 03, x u x 2 , x 3 ) = [Tr( i3 C0 J = i _ o)]^ (E3) 

The [N( 13 CO)/Av] mod and [T R ( 13 CO J = 1 — > 0)] are the column density per velocity 
interval and radiation temperature of the 13 CO J = 1 — > line, respectively, from the model 
fit. Just as there is a family of model curves of r 240 versus Td, the R(a\, a 2 , 03; x\,x 2 , £3) can 
be thought of as a family of curves where each R is given by a particular model r 240 . The 
{cii} are the curve parameters that specify a particular R curve from the family of R curves 
and the {x^ specify, indirectly, a particular point on that curve. Specifically, the {a^} are 

ai = AT = T d -T K , (E4a) 
iV( 13 CO)" 



a 2 = 
and 



Av 



(E4b) 



mod 



«3 = [n(H 2 )] mod . (E4c) 

The {x^ represent a particular observed data point and that point is associated with the 
point on the R versus Td curve that is closest to it, where 

Xl = T d , (E5a) 
x 2 = L,(240/im) , (E5b) 

and 

x 3 = I( 13 CO) . (E5c) 
The uncertainty in the derived column density, cr(y), is given by 

i=i N ' i=i N ' 

The first summation in (E6) represents the contribution of the model-fit uncertainties to 
the computed column density uncertainty. The second summation in (E6) represents the 
contribution of the observational uncertainties to the computed column density uncertainty. 
Specifically, the observational uncertainties result in uncertainties in choosing the nearest 
point on the model curve. To evaluate equation (E6), we will use r = x 2 /x^, where the 
subscript "240" on the r was dropped for brevity. This gives us 

dy OR 

— = — x 3 (E7a) 
acii aa, 
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dy 
dxi 

dy 
dx 2 



dy 

dx 3 



dR 
axi 

dR 

ox 2 
dR dr 
dr dx 2 
dR 
dr 
dR 

tt-xs + R 
dx 3 

OR dr 
dr dx 3 



X3 



x 3 + R 



dR 

= -—r + R 

dr 

Substituting equations (E7a) through (E7d) into (E6) yields 



Ay) 



E 

i=i 



+ 



dR 

dd; 



(dR 
\ dr 



x 3 a (di) 



r 2 - 2rR 



+ 



dR 



dx\ 



x 3 a (xi) + 



dR 

dr 



+ R 2 



Employing 

changes (E8) to 

Ay) = 



a 2 (r)=r 2 [^ + a2 ^ 



x\ 



x\ + 



dR 

dr 



o\x<i) 



(E7b) 



(E7c) 



(E7d) 



(E8) 



(E9) 



Xn 



+ 



R 2 



<J 2 (x 3 ) 



(E10) 



Equation (E10) has the advantage that there is no explicit dependence on <x(x 2 ), dR/dx 2 , 
nor dR/dx 3 . Instead, there is an explicit dependence on r and on cr(r), both of which have 
already been computed for the LTE modeling. Equation (E10) was used in determining the 
uncertainties in the derived column densities. 



We must now compute the derivatives of R. The derivatives with respect to the a,, 
dR/dcii, are determined from the difference in R, AR, that results for a given difference in a 
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specific di, Adi, from the family of model curves. In other words, the following approximation 
is used: 

< E11 > 

dcii Adi 

The uncertainties in the cr(di) were determined by finding the values of a« that resulted in 
the total x 2 rising by a value of xt, since this value was greater than 1 (see discussion of this 
in Section 3.1). Since the model curves were on a grid, changes in Oj occur in discrete jumps 
that increase x 2 by amounts less than or greater than xl- F° r Figure 16, the error bars 
represent the minimum change in Oj that increased x 2 by more than xl- The cr(aj) used for 
deriving cr(y) was the change in a* interpolated linearly between the grid points that gave 
the appropriate increase in x 2 ■ Where there was a limit instead of an upper (lower) error 
bar, an error bar the same size as the lower (upper) error bar was assumed. 

The derivatives dR/dx\ and OR/ Or are with respect to the observed X\ (i.e., Tj) and the 
observed r (i.e., I„(240 /im)/I( 13 CO)). The former derivative is the amount that R changes 
when the observed x\ changes (for fixed r). In other words, when a point in the r versus 
xi plot changes in x 1 , the closest point on the model curve to this point will also change, 
thereby changing R. The latter derivative is similar, but with r in place of x\. Determining 
the closest point on the model curve to the observed data point requires defining distances in 
the r versus xi plane. The distance from any observed point can be defined by normalizing 
changes in x\ to its observational uncertainty, cr(xi), and changes in r to cr(r). Thus we 
define a coordinate system for a given data point as 

X = -^r (E12a) 

Y = . (E12b) 

a(r) 

The distance between the data point and a point on the model curve is simply (X d — X c ) 2 + 
(Yd — Y c ) 2 , where (Xd, Yd) are the data point coordinates and (X c , Y c ) are the coordinates of 
a point on the model curve. Because the curve is a smooth curve, it can be approximated by 
a straight line in an interval that is sufficiently small. The interval of interest is one centered 
on the point on the curve closest to the data point. Varying the X d of the data point, while 
holding Y d constant, or vice versa, then varies the coordinates of nearest point on the curve. 
The variations of these coordinates, (AX C , AY C ), are related to the data point coordinate 
variations by 

AXd = Ax^y? (El3a) 



for constant AYd and 
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^ = A^AVl (E13b) 

for constant AX^. The AX^, AYd are the data point coordinate variations corresponding to 
variations AX C , AY C . Equations (E13a) and (E13b) were derived with simple trigonometry. 
Using equations (E12a) and (E12b), equations (E13a) and (E13b) become 

= (AxOXM + (Ar)^(x,) 
(Axi) c (7 2 (r) 

and 

(Ar)rf " (M)^^) • (E14b) 

The derivatives dR/dxi and dR/dr use the approximations 

ai? Ai? 



dxi (Axi) c 

and 

di? Ai? 



(E15a) 



(E15b) 



9r " (Ar) d 

These derivatives are then computed using equations (E15a) and (E15b) combined with 
(E14a) and (E14b): a particular value of (Axi) c is added to the (x 1 ) c of the point on the 
model curve that is closest to the data point. This implies a particular (Ar) c , which, in 
turn, implies a particular AR. These are placed in equations (E14a) and (E14b) to compute 
the variations induced at the data point, (Axi)d and (Ar)d- These, in turn, are used in 
equations (E15a) and (E15b) to compute the desired derivatives. In the special cases where 
(Axi) c = and (Ar) c = 0, dR/dx± = and dR/dr = 0, respectively, which can be verified 
geometrically and do not come from equations (E14a) and (E14b). 

Generalizing the above to the two-component LVG model is trivial. The definitions of 
the {aj} are now the seven parameters of the two-component models, instead of the three 
parameters shown in equations (E4a), (E4b), and (E4c). 



F. Uncertainty Analysis of the Dust-Derived Column Density for the 

Two-Component Models 

As stated in Section 3.3, the column density of molecular gas derived from the continuum 
observations is derived from expression (40), which is written here as 

Y = xi, x 2 , x 3 ) X(x 1 ,x 2 ) , (Fl) 
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where Y is the observed column density, N d (H), C is the factor multiplying N d i(H), and 
X is Ndi(H). The {xi} are defined as Appendix E, except that x 2 is now T dc . The {a^} 
are the seven parameters that specify each model curve for the two-component models (see 
Section 3.3). The function C is analogous to the R in Appendix E. Consequently, the 
uncertainty derivation is very similar to that in Appendix E. The only other expression 
needed is that for X. This can be in the Wien limit for simplicity, as was done in Section 2.1 
and Appendix D. Specifically, equation (Dl) is now written as 



X(x 1 ,x 2 ) = x 2 D exp — 



(F2) 



As in Appendix E, taking partial derivatives of Y with respect to the {a^} and {x^} is 
necessary. Consequently, we use the following: 



ax = i-f 



and 



dxi 

dX 

dxo 



Xi 



X 

X-2 



(F3a) 



(F3b) 



Taking the partial derivatives of Y with respect to the {a t } and {x^ and adding them as 
in equation (E6) gives 



An = 



1=1 
+ 2 



X 2 + 



dC 



„ dC (7 2 (xi) 

OX i x\ 



dx\ 

dC a 2 (x 2 ) 
dr x\ 



X 2 + C 2 a 2 {X) 



a 2 (r) 



X< 



(F4) 



where r is defined as in Appendix E. The partial derivatives of C were determined using the 
same approach as for the partial derivatives of R in Appendix E. Expression (F4) was used 
to give the uncertainties in the derived N d (H) values for the two-component models. 
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Fig. 1. — The panels above show CO rotational line ratios as a function of kinetic temper- 
ature, T K , in the LTE limit. The solid curves show the ratio of the radiation temperature, 
T R , of the J = 2 — > 1 line to that of the J = 1 — > line. The dashed curves show the 
T R (J = 3 — > 2)/T R (J = 2 — > 1) ratio. The horizontal dotted lines mark the levels of the 
ratios at 80% of their asymptotic values in the high-T K limit. The vertical dotted lines mark 
the T K values for which the ratios achieve this 80% level. Panel (a) shows the behavior of 
these ratios in the case where all the lines are optically thick. Panel (b) shows this behavior 
in the optically thin case. 
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Figure 2a 




Fig. 2.— The DIRBE Surface Brightness Maps for A = 140 and 240 /xm from the 1998-release 
are displayed in the equal-area Mollweide projection with a Galactic longitude- latitude coor- 
dinate grid. The analysis concentrates on the three 10°- 12° circular fields that appear in the 
map as roughly elliptical loops; circles on the sky appear distorted in this projection. The 
fields, from east to west, are the Orion A, Orion B, and A Orionis Fields. Also shown with 
dotted ellipses are the smaller (2° diameter) Orion Nebula and NGC 2024 Fields (which are 
circular on the real sky). The absolute photometric uncertainties of the maps are 10.6% for 
the 140 /im map and 11.6% for the 240 /xm map. All maps have been "cleaned" of zodiacal 
light and had cosecant(\b\) backgrounds subtracted. The contour levels for the 140 /xm map 
are 5, 10, 15, 20, 30, 40, 60, 80, 100, 150, 200, 400, 600, 800, 1000, 1500, lOOOMJysr' 1 and 
for the 240 /tm map are 3, 5, 10, 15, 20, 30, 40, 60, 80, 100, 150, 200, 250, 300, 350, 400, 450 
500, 550, 600, 650, 700, 750 MJy • sr" 1 . 



Fig. 3. — Spatial distribution of the 140/xm/240/xm color temperature, Td, is illustrated 
{(3 = 2 emissivity index unless otherwise stated). Contour levels are 15, 16, 17, 30 K. All 
positions in the Td map where L,(140 pm) < 1 MJy ■ sr _1 or 1^(240 fim) < 1 MJy ■ sr" 1 were 
set to zero, resulting in the close spacing of contours at the edges of the emitting regions. 
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Fig. 4. — Map of the velocity-integrated radiation temperature of the CO J = 1 — > line 
is shown for the Orion region. Contour levels are 0, 0.3, 0.5, 1, 2, 4, 6,...,12K- km-s -1 . The 
thick contour indicates the edge of the mapped area. 
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Fig. 6. — Plots of r 240 or 1^(240 /xm)/I( 13 CO) versus T d are shown for the Orion region. 
The plots include points at more than 5-cr in I„(140/zm), L,(240 /xm), and I( 13 CO). Panel 
a) includes the error bars of the random errors (i.e. systematic errors are not included), 
whereas in Panel b) the error bars are excluded to more clearly reveal the overall pattern in 
the plotted points. The diamonds in this panel are the fiducial points discussed in Section 3.5 
and listed in Table 3. 
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Fig. 7. — The model curves of r 240 versus T d are shown for the LTE case. The plot includes 
the same data points as in Figure 6 with error bars excluded for clarity. Each model curve is 
for a specific value of AT or T K — T d . The curves shown range from AT = —16 K to +16 K 
in steps of 2 K. The diamonds in this figure are the fiducial points discussed in Section 3.5 
and listed in Table 3. 
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Fig. 8. — The model curves of r 240 versus T d that best fit the data are shown for the LTE 
case. The plot includes the same data points as in Figure 6 with error bars excluded for 
clarity. The solid curve represents AT = — 4K and is the best fit to all the data points. 
The dotted curve represents AT = +9 K and is the best fit to only the data points with 
T d > 20 K. The diamonds in this figure are the fiducial points discussed in Section 3.5 and 
listed in Table 3. 




Fig. 9. — The effect of the systematic uncertainties on the resultant AT from the fits of the 
LTE model curves is shown. The effect of these uncertainties was tested by applying the 
scale factors 0.6, 0.8, 1.0, 1.2, and 1.4 to the model curves and fitting the AT for each one. 
The solid line represents the resultant AT values for the fits to all the data (i.e., all the data 
points shown in Figure 6). The dotted line represents the resultant AT values for the fits to 
the data with T d > 20 K. The error bars represent the formal error bars for each model fit. 
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Column Density Comparisons for LTE, One — Component Case 
300 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 




N d (H) (10 20 cm- 2 ) 

Fig. 10. — The plots above compare the gas column densities as derived from gas tracers 
and those derived from dust continuum emission for the case of LTE gas emission and a 
single component. The upper panel shows the gas column density as derived from the 
13 CO J = 1 — > in the LTE, optically thin case versus the gas column density as derived from 
the 140 /im and 240 /im continuum emission. The points in the upper panel plot represent 
the positions where I„(140 /im), 1^(240 /im), and I( 13 CO) are at more than 5-cr, as in Figure 6. 
The lower panel shows the gas column density as derived from the 12 CO J = 1 — > line 
assuming the conversion factor discussed in the text. The points in the lower panel plot 
represent the positions where 1^(140 /xm), 1^(240 /im), and I( 12 CO) are at more than 5-cr. 
Consequently, this sample has more points than in the upper panel: specifically the upper 
panel has 674 points and the lower panel has 1053 points. 
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Fig. 11. — The area filling factors of substructures in the Orion clouds are plotted against 
the molecular gas column densities in the LTE, one-component case. The points in this plot 
represent the positions where 1^(140 //m), 1^(240 /im), and I( 13 CO) are at more than 5-cr and 
where the peak radiation temperature of 12 CO J = 1 — ^ is at more than 3-a. This is a 
total of 372 points. 
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Fig. 12. — The model curves r 240 versus Td for different AT values are shown for the non-LTE, 
one-component case. For the solid curves, the N( 13 CO)/Av value is fixed at 3.2 x 10 15 cm -2 ■ 
(km • s -1 ) -1 and the n(H 2 ) value is fixed at 1 x 10 5 cm -3 — the best fit N( 13 CO)/Av and 
n(H 2 ) values for model fits to all the points in the plot. For the dashed-dotted curves, the 
N( 13 CO)/Av value is fixed at 3.2 x 10 15 cm -2 • (km • s -1 ) -1 and the n(H 2 ) value is fixed at 
5.6 x 10 3 cm~ 3 — the best fit N( 13 CO)/Av and n(H 2 ) values for model fits to only the points 
with T d > 20 K. The curves (solid and dashed-dotted) shown range from AT = — 14 K to 
+ 14K in steps of 2K. The plot includes the same data points as in Figure 6 with error 
bars excluded for clarity. The diamonds in this figure are the fiducial points discussed in 
Section 3.5 and listed in Table 3. 



-85 - 




Fig. 13. — The model curves r 240 versus for different N( 13 CO)/Av values are shown for 
the non-LTE, one-component case. For the solid curves, the AT value is fixed at — 1 K and 
the n(H2) value is fixed at 1 x 10 5 cm" 3 — the best fit AT and n(H2) values for model fits to 
all the points in the plot. For the dashed-dotted curves, the AT value is fixed at —3 K and 
the n(H2) value is fixed at 5.6 x 10 3 cm -3 — the best fit AT and n(H2) values for model fits 
to only the points with Td > 20 K. The solid curves are labelled with their corresponding 
values of N( 13 CO)/Av. The N( 13 CO)/Av values for the dashed-dotted curves follow the 
same order as for the solid curves, i.e., from lower right to upper left. The plot includes the 
same data points as in Figure 6 with error bars excluded for clarity. The diamonds in this 
figure are the fiducial points discussed in Section 3.5 and listed in Table 3. 
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Fig. 14. — The model curves r 240 versus for different n(H 2 ) values are shown for the 
non-LTE, one-component case. For the curves shown, the AT value is fixed at —1 K and the 
N( 13 CO) /Av value is fixed at 3.2 x 10 15 cm" 2 - (km-s" 1 )' 1 — the best fit AT and N( 13 CO) /Av 
values for model fits to all the points in the plot. The best fit AT and N( 13 CO)/Av values 
for model fits to the points with T d > 20 K are similar or identical to those for all the points; 
consequently, those curves would be virtually identical to the solid curves above and are not 
shown. The plot includes the same data points as in Figure 6 with error bars excluded for 
clarity. The diamonds in this figure are the fiducial points discussed in Section 3.5 and listed 
in Table 3. 
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Fig. 15. — The model curves of r 240 versus that best fit the data are shown for the 
non-LTE, one-component case. The plot includes the same data points as in Figure 6 with 
error bars excluded for clarity. The solid curve represents AT = —IK, N( 13 CO)/Av = 
3.2 x 10 15 cm -2 • (km • s -1 ) -1 , and n(H 2 ) = 1 x 10 5 cm~ 3 and is the best fit to all the data 
points. The dotted curve represents AT = -3K, N( 13 CO)/Av = 3.2x 10 15 cm^^km-s -1 )- 1 , 
and n(H 2 ) = 5.6 x 10 3 cm~ 3 and is the best fit to only the data points with T d > 20 K. The 
diamonds in this figure are the fiducial points discussed in Section 3.5 and listed in Table 3. 
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Fig. 16. — The effect of the systematic uncertainties on the resultant parameters from the fits 
of the LVG model curves is shown. The effect of these uncertainties was tested by applying 
the scale factors 0.6, 0.8, 1.0, 1.2, and 1.4 to the model curves and fitting the parameters for 
each scale factor. The left panel shows the resultant AT values, the center panel shows the 
resultant N( 13 CO)/Av values, and the right panel shows the n(H 2 ) values. The solid line in 
each panel represents the resultant parameter values for the fits to all the data (i.e., all the 
data points shown in Figure 6). The dotted line represents the resultant parameter values 
for the fits to the data with T d > 20 K. Notice that the plotted points have been slightly 
displaced horizontally from their true scale factor values for clarity. The error bars represent 
the formal error bars for each model fit and are the minimum grid spacing, for the grid of 
LVG models used, necessary to increase x 2 by a minimum of xl- These formal errors are 
therefore very conservative estimates of the true formal errors. 
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Fig. 17. — The plot above compares the gas column densities as derived from gas tracers 
and those derived from dust continuum emission for the case of non-LTE gas emission and a 
single component. The plot shows the gas column density as derived from 13 C0 J = 1 — > 
using the results from the LVG mehtod versus the gas column density as derived from the 
140 /im and 240 /xm continuum emission. The points in this plot represent the positions 
where 1^(140 /xm), 1^(240 /xm), and I( 13 CO) are at more than 5-cr, as in Figure 6. 
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Fig. 18. — The area filling factors of substructures in the Orion clouds are plotted against 
the molecular gas column densities in the LVG, one-component case. The points in this plot 
represent the positions where 1^(140 /iin), I„(240 /im), and I( 13 C0) are at more than 5-cr and 
where the peak radiation temperature of 12 C0 J = 1 — ^ is at more than 3-a. This is a 
total of 372 points. 
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Fig. 19. — The model curve of r 240 versus T^ c that best fits the data is shown for the non-LTE, 
two-component case. The plot includes the same data points as in Figure 6 with error bars 
excluded for clarity. The parameters and their values are listed in Table 1. The diamonds 
in this figure are the fiducial points discussed in Section 3.5 and listed in Table 3. 
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Fig. 20. — The effect of the systematic uncertainties on the resultant parameters from the 
fits of the two-component, LVG model curves is shown. The effect of these uncertainties 
was tested by applying the scale factors 0.6, 0.8, 1.0, 1.2, and 1.4 to the model curves and 
fitting the parameters for each scale factor. Except for the plots for AT and T d0 , all plots are 
semi-logarithmic where the vertical axes cover about the same logarithmic difference in range 
(about 3 orders of magnitude). This allows easy visual determination of which parameters 
have the smallest systematic uncertainties. 
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Fig. 21. — The plot above compares the gas column densities as derived from gas tracers 
and those derived from dust continuum emission for the case of non-LTE gas emission and 
two components. The plot shows the gas column density as derived from 13 CO J = 1 — > 
using the results from the LVG method versus the gas column density as derived from the 
140 /im and 240 /xm continuum emission. The points in this plot represent the positions 
where 1^(140 /xm), 1^(240 /xm), and I( 13 CO) are at more than 5-cr, as in Figure 6. 
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Fig. 22. — The area filling factors of substructures in the Orion clouds are plotted against the 
molecular gas column densities as computed in the LVG, two-component case. The points 
in this plot represent the positions where I„(140 /jm), 1^(240 /im), and I( 13 CO) are at more 
than 5-0" and where the peak radiation temperature of 12 CO J = 1 — > is at more than 3-er. 
This is a total of 372 points. 
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Fig. 23. — The model curves of r 240 versus that best fit the data are shown for the 
non-LTE, two-component, two-subsample case. The two solid curves represent fits to the 
Tdc < 20 K and Td c > 20 K subsamples. The fit to the former is the curve from a LVG, 
two-component model with c restricted to be > 1. The fit to the latter is an LVG, one 
component model with no such restriction. The dotted curves represent the extensions of 
the model curves across the T dc = 20 K boundary, itself represented as a dashed line. The 
plot includes the same data points as in Figure 6 with error bars excluded for clarity. The 
parameters and their values are listed in Table 2. The diamonds in this figure are the fiducial 
points discussed in Section 3.5 and listed in Table 3. 
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Fig. 24. — The plot above compares the gas column densities as derived from gas tracers 
and those derived from dust continuum emission for the case of non-LTE gas emission, two 
components, and two subsamples. The plot shows the gas column density as derived from 
the 13 CO J = 1 — > using the results from the LVG method versus the gas column density as 
derived from the 140 fim and 240 /im continuum emission. The points in this plot represent 
the positions where 1^(140 /im), L,(240 /iin), and I( 13 CO) are at more than 5-<r, as in Figure 6. 
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Fig. 25. — The area filling factors of substructures in the Orion clouds are plotted against 
the molecular gas column densities as computed in the LVG, two-component, two-subsample 
case. The points in this plot represent the positions where 1^(140 //m), L,(240 /xm), and 
I( 13 CO) are at more than 5-cr and where the peak radiation temperature of 12 CO J = 1 — > 
is at more than 3-cr. This is a total of 372 points. 



- 98 - 




210 " 200 '" 190 



Fig. 26. — Contour maps of the column densities for the LVG, one-component case are shown 
above. The upper panel shows the map for the column densities derived from the 140 /zm 
and 240 /im continuum observations. The contour levels are 3, 5, 7, 10, 20, 30,..., 100, 120, 
140,..., 200 in units of 10 20 H nuclei ■ cm . The lower panel shows the map for the column 
densities derived from the 13 CO J = 1 — > and H I observations. The contour levels are 3, 
5, 7, 10, 20, 30,..., 100, 120, 140,..., 240 in units of 10 20 H nuclei ■ cm" 2 . 
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Fig. 27. — Contour maps of the column densities for the LVG, two-component, two- 
subsample case are shown above. The upper panel shows the map for the column densities 
derived from the 140 /im and 240 jum continuum observations. The lower panel shows the 
map for the column densities derived from the 13 CO J = 1 — > and H I observations. The 
contour levels for both panels are 5, 7, 10, 20, 40, 60, 80, 100, 150, 200,..., 650 in units of 
10 20 H nuclei ■ cm -2 . 
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Fig. 28. — Contour maps of the column densities for the cold (T dl = T K1 < 10 K) dust 
and gas of component 1 in the LVG, two-component, two-subsample case are shown above. 
The upper panel shows the map for the column densities of the cold gas derived from the 
140 /im and 240 /zm continuum observations. The lower panel shows the map for the column 
densities of the cold gas derived from the 13 CO J = 1 — > observations. Along those lines of 
sight with this cold gas the H I column densities were also added. The contour levels for both 
panels are 10, 20, 30,..., 100, 120, 140,..., 200, 300, 400, 500 in units of 10 20 H nucleicm~ 2 . 
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Table 1. Two-Component Model Parameter Values. 



Parameter 


Best-Fit Value 




AT 


OK 




c 


0.04 




Tdo 


18K 




N c0 ( 13 CO) 
Av c 


5 x 10 16 cm -2 • (km • s" 




n c0 


1.8 x 10 4 cm" 3 




Nd( 13 CO) 
Av c 


8 x 10 15 cm -2 • (km • s" 




n c i 


1 x 10 3 cm -3 






5.69 




V 


667 





Note. — The formal uncertainties of AT 
and T^ are both < 3 x 10~ 4 K. The formal 
relative uncertainties of the other parame- 
ters are < 3 x 10" 4 . 
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Table 2. Two- Component, Two-Subsample Model Parameter Values. 

Parameter Best-Fit Parameter Values 

T dc < 20 K Subsample T dc > 20 K Subsample 



AT 



OK 



-3±4K 



Co 
Tao 

N c0 ( 13 CO) 
Av c 

n c o 

N c i( 13 CO) 
Av c 

n c i 



1.0 
18K 

5 x 10 15 cm -2 • (km • s" 1 )" 1 
1.0 x 10 5 cm~ 3 

2 x 10 16 cm' 2 • (km • s" 1 )" 1 (3.2 ± \ 8 ) x 10 15 cm" 2 • (km • s" 1 )" 1 
1 x 10 5 cm" 3 (5.6 ± j "J) x 10 3 cm" 3 



4.60 
525 



9.98 
139 



Note. — The formal uncertainties of AT and T^o for the T^ c < 20 K subsample 
are both < 2 x 10~ 4 K. The formal relative uncertainties of the other parameters 
for this subsample are < 2 x 10~ 4 . 
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Table 3. Inferred Physical Conditions at the Fiducial Points a 



Coordinates 11 


Parameter 15 


Case l d 


Case 2 d 


Parameter Values 
Case 3 e 


Case 4 e 


15.84, 18.26 


AT 


-4 


-1 










N( 13 CO)/Av 


— 


3.2 x 10 15 


(2 x 10 15 ) 8 x 10 15 


(2 x 10 14 ) 5 x 10 15 




n(H 2 ) 


— 


1 x 10 5 


(1.8 x 10 4 ) 1 x 10 3 


(1 x 10 4 ) 5.6 x 10 3 




Distance 


1.17 


1.78 


1.67 


1.92 


17.97, 32.65 


Distance 


0.72 


1.06 


0.31 


0.93 


18.21, 69.46 


Distance 


3.29 


3.43 


0.40 


2.45 


21.25, 11.56 


AT 


+9 


-3 





-3 




N( 13 CO)/Av 




3.2 x 10 15 


(2 x 10 15 ) 8 x 10 15 


3.2 x 10 15 




n(H 2 ) 




5.6 x 10 3 


(1.8 x 10 4 ) 1 x 10 3 


5.6 x 10 3 




Distance 


3.56 


3.67 


2.45 


3.67 


23.35, 48.37 


Distance 


0.99 


0.47 


5.52 


0.47 


26.31, 126.1 


Distance 


4.82 


6.48 


5.29 


5.03 


26.95, 77.87 


Distance 


0.88 


0.45 


1.71 


1.05 



a As depicted in the r 240 versus T d plots of Figures 6, 7, 8, 12, 13, 14, 15, 19, and 23. 



b AT is in units of Kelvins, N( 13 CO)/Av is in units of cm" 2 • (km-s" 1 )" 1 , and n(H 2 ) is in units 
of cm" 3 . "Distance" is the orthogonal distance of the fiducial point from the theoretical curve 
in units of sigmas. The parenthesized quantities in Cases 3 and 4 are the parameter values for 
component 0. In particular, the parenthesized values for N( 13 CO)/Av are really the values for 
the product c * 

c The parameter values for AT, N( 13 CO)/Av, and n(H 2 ) are not explicitly listed for some 
fiducial points. The parameter values for these points are the same as those of the points with 
the listed values [e.g. point (15.48, 18.26) or (21.25, 11.56)] that lie above the points without 
the listed values. 

d Case 1: LTE, one-component models. Case 2: LVG, one-component models. 

e Case 3: LVG, two-component models. Case 4: LVG, two-component, two-subsample models. 
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Table 4. Gas Masses Sampled by 13 CO, 12 CO, & HI (M ) 



Field a 


Case l b 


Case 2 C 


Case 3 d 


Case 4 e 


12 CO & H I f 


Orion Nebula 


21000 


18000 


35000 


31000 


20000 


Orion A 


88000 


85000 


146000 


129000 


122000 


NGC2024 


10000 


9100 


18000 


9200 


17000 


Orion B 


93000 


88000 


248000 


189000 


104000 


A Orionis 


38000 


36000 


92000 


57000 


45000 


TOTALS 


220000 


209000 


486000 


376000 


271000 



a See W96 for definitions of the fields. 

b Using LTE, one-component models and observations of 13 CO and 
HI. 

c Using LVG, one-component models and observations of 13 CO and 
HI. 

d Using LVG, two-component models and observations of 13 CO and 
HI. 

c Using LVG, two-component, two-subsample models and observa- 
tions of 13 CO and H I. 

f Adopting an X-factor of 2.6 x 10 20 cm~ 2 • (K ■ km ■ s -1 ) -1 . 

g Sum of Orion A, Orion B, A Orionis fields. 
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Table 5. Gas Masses Inferred from 140, 240 /zm Continuum Emission (M ) 



Field a 


Cases 1 & 2 b 


Case 3 C 


Case 4 d 


Orion Nebula 


14000 


29000 


27000 


Orion A 


107000 


169000 


151000 


NGC2024 


9700 


14000 


9700 


Orion B 


81000 


225000 


179000 


A Orionis 


63000 


122000 


82000 


TOTAL 6 


251000 


516000 


411000 



a See W96 for definitions of the fields. 

b Using one-component models and continuum 
observations. 

c Using two-component models and continuum 
observations. 

d Using two-component, two-subsample models 
and continuum observations. 

c Sum of Orion A, Orion B, A Orionis fields. 
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Table 6. Best Estimates of Total a Gas Masses (M ) 



Tracer 


One-Component Case b 


Two-Component Case c 


Derived from Continuum 


251000 


411000 


Derived from 13 CO J = 1 — » 


d 263000 


429000 


a The total hydrogen mass 


from the 3 big fields (i.e., 


sum of Orion A, Orion B, 



A Orionis fields. 

b Using LVG, one-component models. 

c Using LVG, two-component, two-subsample models. 
d After correcting for insufficient coverage. 
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Table 7. Masses a of Gas in a Cold (T K < 10 K) Component (M ) 



Field b 


Case 3 C 


Case 4 d 


Orion Nebula 


14000 


14000 


Orion A 


43000 


28000 


NGC2024 


5700 


6600 


Orion B 


155000 


126000 


A Orionis 


54000 


28000 


TOTAL c 


252000 


181000 



inferred from 13 CO and HI ob- 
servations. 

b See W96 for definitions of the 
fields. 

c Using LVG, two-component 
models and observations of 13 CO 
and H I. 

d Using LVG, two-component, 
two-subsample models and obser- 
vations of 13 CO and H I. 

e Sum of Orion A, Orion B, 
A Orionis fields. 



